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ABSTRACT 


An algorithm is developed and evaiuated for discriminating among clouds, snow 
cover and clear land. The multispectral technique uses daytime images of AVHRR 
channels 1 (0.63pm). 3 (3.7m) and 4 (11.0fm). Reflectance is derived for channel 3 
bv using the channel 4 emission temperature to estimate and remove the channel 3 
thermal emission. Separation of ciouds from snow and land is based primarily on this 
derived channel 3 reflectance. Using this technique, observed reflectance in channel 3 
is 2 to 4 percent for snow, 3 to 10 percent for land, 2 to 27 percent for ice clouds and 8 
to 56 percent for liquid clouds. These values overlap for thin cirrus and snow, so the 
routine then attempts analysis of cirrus based on its different transmissive properties 
between channels 3 and 4. Six images were analyzed and the total cloud cover was 
verified against a total of 110 conventional surface observations using the standard 
categories of clear, scattered, broken and overcast. The routine was quite successful, 
with the analyzed sky cover being within category for 55 percent of the stations, one 
category different for 33 percent, 2 categories different for 9 percent and 3 categories 
different for 3 percent of the stations. A major remaining problem is discrimination 
between ice clouds and snow cover due to the great similarity of reflective properties of 


these two surfaces. 
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10 Nov 86 subscene 2 - NOAA-9 AVHRR channel |, full resolution 
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I. AUTOMATED CLOUD ANALYSIS 


A. INTRODUCTION 

One of the critical factors in the conduct of military operations is knowledge of 
the environment in which they are to take place, particularly regarding cloud cover. 
Clouds can dramatically affect such missions as aerial intelligence gathering, air 
refueling and tactical employment of weapons that use the visible-infrared portion of 
the electromagnetic spectrum to designate targets. Cloud cover analyses can be 
accomplished with synoptic weather observations but data-void areas severely limit 
their usefulness. Techniques have been developed to produce cloud analyses using 
satellite imagery, where human image analysts make subjective interpretations of cloud 
cover based on brightness and texture contrasts. Manual analysis takes a great deal of 
time to accomplish and quick response for large geographical areas in critical situations 
is very difficult. The Air Force Global Weather Central (AFGWC) developed an 
automated cloud analysis system in order to introduce objectivity to the analvsis 
procedure based on sound physical principles. and also to address the Air Force's 
mission requirement for rapid cloud analyses. The 3-Dimensional Nephanalysis 
(3D NEPH) svstem became operational in January 1970. This pioneering effort did not 
have the benefit of experience from other systems. Various improvements have been 
made over the years, and in 1984 a second-generation model called the Real-Time 
Nephanalysis (RTNEPH) became operational. It is the only known automated system 
capable of processing and interpreting the tremendous volume of satellite data that are 
available today and integrating them with conventional information to provide a high- 
resolution. three-dimensional cloud analysis data base for the entire world. 

In addition to global analysis models, efforts have been made to develop 
automated analysis routines for minicomputers that can be used in a theater or 
battlefield environment. Wash er al. (1985) developed a cloud and precipitation 
analysis program for an interactive minicomputer system which uses geostationary 
infrared and visual data. This type of analysis package could be used in a weather 
station or a tactical van to provide cloud cover analyses and forecasts to operational 


commanders on a real-time basis. 


Ih. 


B. SNOW/CLOUD DISCRIMINATION 

Various long-term problems have plagued automated systems, one of the most 
significant being discriminating low clouds from snow cover. Middle and high clouds 
can be detected against snow cover due to the thermal contrast detectable at infrared 
wavelengths, but in the regions of the electromagnetic spectrum used most often in 
satellite analysis, low clouds and snow cover have similar radiometric properties. In 
the visible portion of the spectrum both have high albedos, and in the infrared portion 
both have similar thermal properties. Thus there is little contrast in images at these 
spectral wavelengths. 

Bunting ef al. (1977) discussed the reflectance properties of snow and clouds at 
visible and near-infrared wavelengths, using data from the Earth Resources Experiment 
Package which flew aboard Skylab. They analyzed the imagery in pairs: one set in the 
visible spectrum and the other in the near-infrared spectrum (1.554m to 1.75u4m). The 
two sets showed high reflectance for snow, water clouds and ice clouds in the visible 
spectrum, but in the near-infrared the reflectance was high for water clouds, medium 
for ice clouds and very low for snow. They concluded that these differences could be 
the basis for an automated routine which makes an analysis decision based on the ratio 
of visible reflectance to near-infrared reflectance. Bunting and d’Entremont (1982) 
tested six automated classifiers on data from a special Defense Meteorological Satellite 
Program (DMSP) sensor in the wavelength band from 1.51pm to 1.634m, along with 
visual and infrared data. They used the classifiers successfully to distinguish water 
clouds, ice clouds, snow cover and other cloud-free surfaces. Their technique was so 
successful that plans are now being made to add this sensor as an operational 
component of the DMSP satellite. 

Kidder and Wu (1984) showed that there is a contrast in brightness temperature 
between snow and low clouds in daytime imagery of channel 3 (3.7m) of the National 
Oceanic and Atmospheric Administration (NOAA) polar-orbiting Advanced Very High 
Resolution Radiometer (AVHRR). They related this contrast to solar reflection, which 
is a substantial component of the channel 3 radiance measurement along with thermal 
emission. Since low clouds and snow have similar emission temperatures, this contrast 
is the result of a difference in their solar reflection at 3.7m which, in theory, is just a 
few percent for snow and about 20 percent for clouds. This results in a warmer 


brightness temperature for clouds (darker image) than for snow (whiter image). 


[2 


Fye (1978) documented the 3DNEPH svstem and discussed the problem of 
snow,;cloud discrimination. Cloud cover is overestimated in areas of new snow and in 
areas where the ice analysis contains too little ice. Underestimation or no analysis of 
clouds results when snow melts rapidly or the ice analysis contains too much ice. 
Presently, the system uses the AFGWC Snow Cover Model and grid points that have 
snow or ice as a background are treated as missing for the cloud analysis routine. A 
solution to this problem has not yet been implemented on the RTNEPH system. This 
has a significant effect on the analysis in the late fall through early spring because 
about 30 to 35 percent of all the model grid points in the northern henusphere are 


snow and/or ice covered. 


Can UitestS OBJECTIVES 

The first objective of this thesis is to calculate and report solar reflectance in 
AVHRR channel 3 using a method that estimates the thermal emission from channel 4 
and removes this portion from the channel 3 radiance measurement. Derived values 
will be compared to theoretical 3.7m reflectance for snow cover, land and clouds. 
During the 1990’s both NOAA and the Department of Defense plan to launch 
satellites with a 1.6f1m sensor for daytime snow/cloud discrimination, and the NOAA 
satellite will have a 3.711m sensor for use at night. Even with this change, there are 
good reasons to pursue this work now. First, by the time the new sensor flies there 
will be over 10 years of archived daytime 3.7jtm imagery. Second, since snow and 
cloud reflectances are similar at 1.641m and 3.7m, a data set of reflectances from 
channel 3 would be very useful in testing software for the new sensor. 

The second objective of this thesis is to use calculated channel 3 reflectance to 
develop an automated cloud analysis routine that separates snow cover, clouds and 
snow-free land. It will be designed to handle various combinations of cloud and 
surface features, but the focus is on the separation of low clouds and snow. The 
routine will not rely on sophisticated statistical analysis or on artificial image 
enhancement techniques, but rather on basic physical principles of radiative transfer in 
the earth’s atmosphere. In addition to simular visible reflective properties, low clouds 
and snow have nearly identical infrared thermal properties. Because of this and the 
results of the previously mentioned studies, the routine 1s based primarily on the 
reflective properties of these surfaces in AVHRR channel 3. A multispectral technique 
is developed and evaluated using visible and infrared data along with channel 3 to 


process davtime AVHRR images. 
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In Chapter II the necessary theoretical background 1s presented. The details of 
satellite data processing are presented along with a theoretical discussion of the 
reflection and enussion characteristics of the cloud and surface features to be analyzed. 
A discussion of the satellite data and the techniques of data analvsis are presented in 
Chapter Ii] along with observed values of channel | and channel 3 reflectance for 
these surfaces. Finally, the analysis routine is developed using these measured 
reflectance values. Chapter IV discusses the performance of the analysis routine when 
applied to six images containing various combinations of snow and cloud cover, and 
sun-satellite geometry. Data were taken from two different NOAA satellites to 
illustrate the general utility of the routine and to provide additional statistical 
independence to the verification of the analyses. Total sky cover on the analyses is 
verified against the available surface observations using the standard conditions of 
clear, scattered, broken and overcast. The reported statistics include the percentage of 
stations where the analyzed sky cover matched the observed sky cover, and the 
percentage of stations with a one-category, a two-category and a three-category 


difference between the analyzed and observed sky cover. 
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Il. THEORETICAL FOUNDATION 


A. DISCUSSION 

In general terms, detection of cloud cover on a satellite image depends on the 
contrast between the clouds and their background. At visible wavelengths this contrast 
is in terms of reflectance differences, and at infrared wavelengths this contrast is in 
terms of differences in brightness temperature. At middle infrared wavelengths 
(AVHRR channel 3) the contrast is a mixture of difference in reflectance and difference 
in brightness temperature. This chapter presents the theoretical foundation for the 
algorithm developed in this thesis. Included here is a theoretical discussion of the 
reflective and thermal characteristics of the features to be analyzed, the satellite data 
processing to derive these properties and a preliminary assessment concerning how the 
routine might work based on the theoretical characteristics of the features. See 
Appendix A for a list of symbols and constants used. 

Two visible and two thermal infrared channels are available on the AVHRR 
instrument. In addition to these, channel 3 is at a wavelength (3.7m) that has 
contributions from both thermal emission and solar reflection. Channel 1 (0.63ftm) is 
completely in the visible spectrum while channel 2 (0.87f1m) extends into the near- 
infrared portion of the spectrum. Channel | is the visible channel used in this study. It 
will be shown later that the method of deriving the reflectance in channel 3 depends on 
channel 4 (11.0,m), so it is the infrared channel used in this study. 

Throughout this thesis, the satellite-measured radiance during daytime 1s 


approximated as follows: 
eet)! + 1(0.0.@)lcosd al) 


The first term on the right hand side is the contribution to the measured radiance from 
thermal emission of the viewed surface assuming the transmissivity between it and the 
satellite is 1.0. The amount of radiance from thermal emission reaching the satellite 1s 
determined by the emissivity (€) of the viewed surface. The Planck function relates the 
emitted monochromatic intensity with the wavenumber and temperature of the 


emitting surface as follows: 


2he?v> 
(  — ————— (2) 
exp(hcv/KT) - | 


The second term on the right hand side of Eq. 2.1 1s the contribution to the measured 
radiance due to solar reflection. The amount of reflected solar radiance reaching the 
satellite is determined by the incident solar radiance (1) which is weighted by the cosine 
of the solar zenith angle (.). and the reflectance (r}) of the viewed surface. The 
directional reflectance 1s a function of the solar zenith angle, the satellite zenith angle 
(9) and the horizontal angle between them (@). This sun-satellite geometry is 
illustrated in Fig. 2.1. Isotropic reflectance is related to the directional reflectance by 
the anisotropic reflectance factor (f). This relationship 1s expressed as (Taylor and 
Stowe, 1984): 


_ 1(8,,8.9) 
a ae i‘ 
The anisotropic reflectance factor arises due to the directional dependence of the 
scattering mechanisms. It is the ratio between the radiant exitance in a given direction 
assuming the surface reflects isotropically, and the actual radiant exitance. For 
example, a radiance measurement in a given direction with a factor of 1.0 gives the 
correct radiant exitance, while a factor of 1.5 with the same measurement means the 
isotropic assumption would yield a 50 percent overestimate of radiant exitance. 

It is assumed that any viewed clouds are optically thick so that upwelling 
radiance from below them does not reach the satellite. This assumption breaks down 


when considering thin cirrus clouds and this will be discussed later. 


B. REFLECTANCE 
This section reviews the results of various theoretical studies concerning the 
reflective properties of the surfaces to be analyzed. These values then will be used for 
comparison to observed values in the data analysis section of Chapter ITI. 
1. Liquid Clouds 
The reflectance of liquid clouds is dependent upon cloud optical thickness and 
sun-satellite geometry. Tables | and 2 show the results of reflectance calculations from 
the model of Shettle and Weinman (1970), which is based on Mie theory and the delta- 
Eddington approximation. The values show the expected extremes of reflectance for 


clouds of various physical properties which are coded in the tables as follows: 
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Fig. 2.1 Earth-sun-satellite geometry. 


1. First digit = liquid water content (x 107! g/m?). 

2. Second digit = droplet mode radius ( Xx 10: yim). 
For example, a code of 66 means the liquid water content 1s 0.6g/m° and the droplet 
mode radius is 0.6j1m. The extinction coefficient, asymmetry factor and the single 
scattering albedo were used as input to the computer program along with the cloud 
optical depth and solar zenith angle to calculate the reflectance. 

Table | shows the strong dependence of reflectance at visible wavelengths on 
cloud thickness. To a lesser degree this reflectance is dependent on liquid water 
content and solar zenith angle. The same gencral characteristics for 3.7J1m are 
apparent in Table 2, but in this case the dependence on liquid water content and cloud 
thickness is not as great as the dependence on solar zenith angle. Absorption by liquid 
Water at 3.7ftm is so great that a cloud has essentially an infinite optical depth with a 
thickness of just a few tens of meters. These values agree closely with those reported 
by Hansen and Pollack (1970), and Welch ez al. (1980). The major differences in the 


values listed in Table 2.6 of Welch er al. and the values listed here in Tables | and 2 


ie 


TAB VED 


Calculated reflectance of water clouds at 0.63.1m. 


Solar Zenith 
Cloud 
Code Thickness 45 5S 65 is 85 
34 100m 0; 395 0.43] 0.500 0.577 0.656 
1000m 0.86] 0.875 0.890 0.907 0.924 
36 100m 0.269 0.327 0.404 0.499 0.596 
1000m O,794 0.814 0.836 0.86] 0.887 
ai 100m 0.209 07262 0.340 0,446 0.561 
1000m O35 0.760 0.790 0.82] 0.854 
64 100m 0.549 0.593 0.642 0.696 0.753 
1000m 0.926 0.933 0.94] 0.950 0.959 
66 100m 0.428 0.482 0.545 0.614 0.686 
1000m 0.885 0.897 0.909 0.923 0.937 
68 100m 0.349 0.406 0.478 0.559 0.641 
1000m 0.848 0.863 0.879 0.897 0.916 
94 100m 0.648 0.682 0.72] 0.763 0.807 
1000m 0.949 0.954 0.960 0.966 0.972 
96 100m Osa US) 0.628 0.684 0.743 
{000m 0.92] 0.928 0.937 0.947 0.956 
98 100m 0.448 0.500 0.56] 0.628 0.697 
1000m 0.893 0.904 0.915 0.928 0.94] 





are due to the different solar zenith angles and a slight difference in wavelength. These 
calculated values, along with the measured values which will be reported in Chapter 
III, are used to determine the thresholds for separating hquid cloud surfaces from other 
surfaces. 
2. Ice Clouds 

The reflectance of ice clouds is not only dependent on the optical thickness 
and sun-satellite geometry as with liquid clouds, but also on complicating factors of ice 
particle shape and size distribution. Modeling the reflectance of ice crystals is an 
extremely difficult task because of these complications. Welch et al, (1980) discussed 
the physical characteristics and properties of ice crystals. They reported reflectance at 
3.341m, which is slightly outside the channel 3 spectral band, for a wide range of shape 


and particle size distributions. They indicated a range of values of reflectance at 3.3p1m 
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Weis, 2 


Calculated reflectance of water clouds at 3.7m. 


Solar Zenith 


Cloud 
Code Thickness 45 55 65 75 85 
34 100m 0.244 0.280 0.326 0.384 0.456 
1000m 0.250 O55 0.330 Uo 0.459 
| 36 100m 0.159 0.192 0.235 270 0.361 
| 1000m 0.170 0.201 O42 0.295 0.365 
38 100m OalhZ 0.14] 0.180 Ci 0.298 
1000m (al 25 0.152 0.189 0.237 0.302 
64 100m 0.250 0.285 0.330 0.387 0.459 
1000m 0.250 0.285 0.330 0.387 0.459 
66 100m 0.169 0.200 0.241 0.295 0.364 
1000m 0.170 0.201 0.242 0.295 OCS 
| 68 100m Orin Qabol 0.188 0.237 0.302 
1000m Oris Oe 0.189 O237 0.302 
94 100m 0.250 0.285 0.330 0.387 0.459 
1000m 0.230 0.285 0.330 0.387 0.459 
96 100m 0.170 0.201 242 0.295 O05 
1000m 0.170 0.201 O47 0.295 0.365 
9§ 100m Oei2> On Se 0.189 O.257 0.302 
1000m Ont 25 Oe sive 0.189 e238 / 0.302 





of 0.001 to 0.007 with a solar zenith angle of 0°. They also calculated reflectance at 
O.55f1m and 0.764tm. The values at both of these wavelengths are similar and can be 
used as an estimate of reflectance in channel | which 1s between these two wavelengths. 
Thev reported a range of reflectance values for 0.554m and 0.76m between 0.01 and 
0.93 with most values in the range from 0.60 to 0.90. Arking and Childs (1985) 
reported channel 3 reflectance of ice clouds as a function of optical thickness between 
0.1 and 100.0, and of ice particle size between 2m and 32pm. Their values of channel 
3 reflectance range between 0.01 and 0.30 for a solar zenith angle of 60°, a satellite 
zenith angle of 41.2°, and a relative azimuth angle of 110°. Finally, Bunting (personal 
communication, 1986) reported an average channel 3 reflectance of 0.05 for ice clouds. 
Reflectance of ice clouds will be measured and reported in Chapter III, and the 
measured and theoretical values will be used to establish thresholds for separating ice 


clouds from the other surfaces. 


ie 


3. Clear Land 


The reflectance of clear land is dependent upon soil type, vegetation type and 


coverage, and moisture content. Bunting and Hardy (1984) reported the albedo of 


snow-free land at visible wavelengths to be approximately 0.10. Bunting (personal! 


communication, 1986) reported an average value of 0.15 at 3.71tm. 


4. Snew Cover 


The reflectance of snow is dependent upon a variety of factors: snow grain 


size, liquid water content, solar zenith angle, cloud cover, snowpack thickness and 


impurities. Wiscombe and Warren (1980) developed a model for the spectral albedo of 


snow and considered these factors: 


ic 


bo 


Snow grain size or age. The albedo of snow decreases as it ages due to an 
increase in absorption and more forward scattering of larger grains which 
develop from the aging process. At visible wavelengths, the reduction is only 
about 10-15 percent, but at 3.71t1m this reduction is a factor of 2 or more for 
grain radi larger than 100}1m. 


Liquid water content. Snow melt affects albedo because water replaces air 
between individual snow grains. This has the effect of increasing the effective 
snow grain size. The reduction in albedo is not as great as the effect due to 
aging. 


Solar zenith angle. Albedo increases at all wavelengths as the zenith angle 
increases. Relative to an overhead sun, the albedo increases only by a few 
percent at visible wavelengths, but at 3.7m the albedo increases from less than 
0.005 to 0.025. These calculations are for a flat snow surface. When the snow 
pack 1s uneven due to hills and valleys, the solar zenith angle is a function of 
location, causing the albedo to vary within a small area. This may occur on a 
large enough scale to appear as a variation of brightness from pixel to pixel on 
a satellite image. 


Effects of cloud cover. Cloud cover influences spectral albedo by changing 
direct radiation into diffuse radiation. This should be important only when 
there is thin cloud cover through which a satellite can see snow cover. 


Snowpack thickness. If the snow cover is thin, the underlying surface affects 
the albedo. This effect is wavelength dependent and also dependent on the 
snow grain size. Because of the strong absorption at 3.7m, less than 1 mm of 
depth is required to make the snow cover optically semi-infinite. At visible 
wavelengths the radiation can penetrate to greater depths and in general, a 
depth of 20 cm is required in order not to “see” the ground. 


Wiscombe and Warren (1980) used Mie theory for single scattering and the delta- 


Eddington approximation for multiple scattering. Their model allows snow albedo 


calculations for the entire solar spectrum (0.3j}tm SA 5.0m) for diffusely or directly 


incident radiation at any zenith angle. Calculations of snow albedo were made using 
this model and compared with observations, and there is good agreement at 3.7ptm. 
Values used for comparison with observations in this study are 0.005-0.025. This range 
of values is applicable for a semi-infinite snow pack with a grain radius of 504m to 
200um, and a range of solar zenith angles between 40° and 80°. At. visible 
wavelengths, there is a considerable difference between calculated and measured values 
wnich thev attributed to impurities in the snow. This is an important consideration for 
this study since the cases are taken sufficiently close to the industrialized regions of the 
United States and Canada. The value for 0.63um should be 0.75 to 0.80 for semi- 
infinite snow cover with a solar zenith angle of 60°. This value is less for a smaller 
zenith angle or snow cover of less than 20 cm in depth, and higher for a greater solar 


zenitn angle. 


C. INFRARED TRANSMISSIVE PROPERTIES 

As shown above, snow and optically thin ice clouds have similar reflective 
properties, thus another method must be used to separate them. It is expected that, 
since the channel 3 reflectance of optically thick ice clouds is greater than that of snow, 
the thresholds in the analysis routine can be tuned finely enough to correctly analyze 
most of the ice clouds, but there will be cases for which this fails. The easiest 
approach would be to analyze as much of the data as possible in an effort to find a 
threshold emission temperature between snow and ice clouds. However, this does not 
allow for the fact that snow is warmer in the fall and spring than in the middle of 
Winter. and also for the fact that snow temperature 1s a function of latitude. In reality, 
the snow temperature is too variable to allow an analysis routine to use a single 
threshold value. One could get around this problem if the surface temperature for each 
pixel were known, but the analysis routine will be designed to work with the satellite 
data alone so this is not possible. 

Another possibility is to take advantage of the difference in the transmissive 
properties of ice clouds between channel 3 and channel 4. The assumption made at the 
beginning that no radiance from below the clouds reaches the satellite breaks down in 
the case of thin cirrus. However, this can be exploited to differentiate between cirrus 
and snow because of a difference in transmissivity between channels 3 and 4. Eg. 2.1 
for the channels 3 and 4 radiance measurement becomes: 


L, = 1,€,(s)B,(T,) + €3(c)B,(T,) + 1,(6,,0.9)Icos8, (2.4) 


3 


21 


and 
L, = T4£,(s)B,(T.) oP E,(c)B,(T.), (25) 


where T is the transmissivitv of the cirrus cloud. 

Neglecting solar reflection and assuming the clouds are optically thick, the 
measured radiance for both channels includes only the thermal emission from the 
cloud, and inverting Eq. 2.2 to get a brightness temperature yields the same result for 
both channels. In the case of thin cirrus, and again neglecting solar reflection, there is 
a cifference in its transmissive properties between channels 3 and 4, which results in a 
different amount of surface radiance for each channel getting through the clouds to be 
sensed by the satellite. Stephens (1981) argued that, first, the Planck blackbody 
emission as a function of temperature varies differently for these two channels and, 
second, cirrus transmits more surface radiation at 3.7m than at 1lfm. When cirrus is 
present, the channel 3 cloud emission is small compared to surface emission, making 
the satellite-measured radiance for this channel primarily due to surface radiation 
transmitted through the cirrus cloud. In channel 4 there is significant emission from 
within the cloud compared to surface emission and the cloud transmits less of the 
surface radiation. The satellite-measured radiance in channel 4 is then primarily due to 
emissicn from the cloud. Both of these effects combine fapproduce a dilleseaeeum 
brightness temperature. 

In reality, the effects of solar reflection cannot be neglected. Solar reflection is 
larger in channel 3 than in channel 4 and this contributes to the difference in brightness 
temperature. The magnitude of solar reflection for cirrus is comparable to the 
magnitude of thermal emission, thus solar reflection contributes a great deal to the 
difference in brightness temperature. However, there are certain conditions under 
which the reflectance of snow cover and ice clouds is equal, so when comparing Eq. 2.4 
for these surfaces, the third term is the same for both. Under these conditions the 
distinction between them arises because T, is not equal to Tt, for cirrus, but T, 1s very 
nearly equal to t, for the clear sky case. Fig. 2.2 shows this difference in brightness 
temperature for cirrus clouds from model calculations of Stephens (1981). The dotted 
line 1s a nighttime case, which shows the effects of thermal emission alone, and the 
solid line is a daytime case with an overhead sun. In this study the solar zenith angle 


is greater than 45° so the temperature difference due to reflection is less, but still 


significant. At very low optical depths the temperature difference is nearlv all due to 
solar reflection and the upwelling radiance is essentially all from the surface. The 
temperature difference increases as optical depth increases and, for the nighttime case, 
reaches maximum at an optical depth of about 2.0. At higher optical depths the 
temperature difference continues to increase, but the effect of transmissivity differences 
decreases to zero. At optical depths higher than 10.0 the difference in brightness 
temperature is due only to solar reflection. 

The distinction between cirrus and snow cover (clear skv) can be made only when 
the effect of the difference in its transmissive properties between channels 3 and 4 
outweighs the effect of solar reflection. For example, assume a cirrus cloud is to be 
analyzed that has a channel 5 reflectance of 0.03. There is a total difference in 
brightness temperature of 10K, half due to reflection and half due to to transmissive 
differences between channels 3 and 4. Snow cover with the same reflectance produces 
a temperature difference of 5K so in this case it is distinguished from cirrus. If the 
snow reflectance is higher than the cloud reflectance, the temperature difference is also 
higher and this would eliminate some or ail of the distinction present when the snow 
reflectance is the same as that of cirrus. If the effect of transmissivity differences had 
contributed 8K to the total temperature difference in this example, only 2K would be 
subject to being eliminated by variability of solar reflection. As seen in Fig 2.2, the 
effect of transmissivity differences is larger than the effect of solar reflection only for 
optical depths between about 0.5 and 2.0. It remains to be seen how many cirrus 
cloud pixels will be handled correctly based on channel 3 reflectance alone and thus 
how many of the remaining cirrus clouds will fall within this optical depth range. 

For this study, a unitless factor is derived from this difference in brightness 
temperature for use by the analysis routine. Development of this temperature factor 
begins with normalizing the temperature difference by the channel 4 brightness 


temperature: 


Hales 


2.6 
= (2.6) 


= 
4 


It follows that, using this approach, the temperature factor is proportional to the ratio 


of the channels 3 and 4 brightness temperatures: 
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Fig. 2.2 Difference in channel 3 and channel 4 brightness temperature 
as a function of cloud optical depth. Solid line is daytime case with overhead sun, 
dotted line is nighttime case (adapted from Stephens, 1981). 
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These numbers range from approximately zero for the clear sky case to values between 
zero and 1.0 when thin cirrus is present. It is nearly impossible to choose a threshold 
from this small range of values, so the reciprocal of Eq. 2.7 was used to choose a 
threshold. This means that, in qualitative terms, for the clear skv case the temperature 
difference is small, making the factor large. For thin cirrus the temperature difference 


is larger, making the factor smaller than with the clear skv case. 


D. SATELLITE DATA PROCESSING 
This section discusses how the satellite data were processed in order to derive the 


reflective and thermal properties presented above. 


24 


1. Channel 1 

AVHRR Channel | data are preflight calibrated in terms of albedo (Lauritson 
et al., 1979). The measured directional reflectance was converted to isotropic 
reflectance by weighting the measurement by the cosine of the solar zenith angle and 
the anisotropic reflectance factor as discussed above. This reflectance is very similar 
for snow cover and clouds, but it can be used to separate land from clouds and snow 
EOucr, 

2. Channel 3 

Channel 3 AVHRR data are calibrated in units of mW/m?-srecm’!, and during 
Caytime the radiance measurement in this channel is comprised of both solar reflection 
and thermal emission. Fig. 2.3 illustrates the relative contributions of each using 
incident solar irradiance measurements of Thekaekara er al. (1969), assuming a typical 
cloud reflectance of 0.05 to 0.40 and assuming a Planck blackbody emission 
temperature of 260K. The dashed vertical lines in the figure define the approximate 
location of the 50 percent response function boundaries of NOAA-9 channel 3. Fig. 
2.3 shows that solar reflection is the greatest contributor for most clouds, while for the 
lowest reflecting surfaces of ice clouds, land and snow cover the two contributions are 
nearly equal. Bunting (1986, personal communication) suggested the method discussed 
here for removing the thermal emission in channel 3. There are two major assumptions 
invelved in using this technique (Ruff and Gruber, 1983). First, atmospheric 
transmission effects are neglected. The most important absorbing constituent in the 
atmosphere in this portion of the spectrum 1s water vapor and in this study the data 
were taken from high latitude, cold air masses where there its little water vapor present. 
Second, it was assumed that the emissivitv of a viewed scene is 1.0. Thermal emission 
for the scene was then determined directly from the channel 4 radiance measurement. 
Error results from this assumption due to the fact that the emissivity of any surface 1s 
variable between channels 3 and 4 and not equal to 1.0. It 1s impractical to try to 
determine the emissive properties of the viewed scene so this effect was not corrected. 
The error resulting from these assumptions should have the greatest effect when 
deriving the reflectance for snow cover, land, or ice clouds because the magnitude of 
thermal emission for these surfaces is comparable to that of solar reflection. The least 
effect from the resultant error should occur with hquid clouds since for these surfaces 
the magnitude of solar reflection is much greater than that of thermal emission (see 
pice 5). 
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Fig. 2.3 Refiected solar intensity and terrestrial thermal intensitv for channel 3. 
Reflectance is 5 and 40 percent, solar zenith is 60° (solid lines), 
blackbodv temperature is 260K (dotted line). 


Pon channel. 24, 2-1 becomes: 
ea JI) 8 r,(8.8,9)I,cos0. | (2333) 


See Appendix B for the calculation of the channel 3 incident solar radiance. Thermal 
emission for channel 3 was estimated by solving Eq. 2.1 for the Planck function based 


on the cnannel 4 measured radiance: 
BAT) = Lit, (29) 


The temperature was found bv inverting the Planck function, and was then used in Eq. 


2.2 to find the thermal emission in channel 3. 


Ignoring atmospheric absorption, emittance for channel 3 is related to 


reflectance by: 


eros, (2.10) 
In Eq. 2.10 isotropic radiation was assumed. Substituting Eqs. 2.3 and 2.10 into Eq. 


2.Segives: 
ea) t,)B,(1) r,(1,cos®. )f (2 lel) 
Finally. solving for reflectance: 


3 (i,cos8, )f- B,(T) ee) 


Eq. 2.12 was used to derive the channel 3 reflectance and this calculation is the central 
component of the analysis routine. 

To calculate the brightness temperature, the data were calibrated and then Eq. 
2.2 was inverted to solve for temperature. The central wavenumber used in this 
calculation is based on a range of temperature values between 225K and 320K 
(Lauritson er al., 1979), which should be a valid range of temperatures for the time of 
vear from which the data were taken. 

3. Channel 4 

Channel 4 AVHRR data are calibrated in units of mW/m?.srecm’!. A non- 
linear correction (Lauritson et al., 1979) must be added when converting from radiance 
to temperature. A polynomial was fit to the published corrections for NOAA-9 for use 
in the computer program (Gerald and Wheatley, 1984). The polynomial derived for the 
NOAA-10 corrections was unstable so linear interpolation was used instead. For ease 
in computer processing all channel 4 temperatures were assumed to be in the range 
from 225K to 275K so that only one central wavenumber need be used in the 
calculations (Lauritson er al., 1979). This range of temperatures provides adequate 
separation of the various features to be analyzed. Temperatures less than or greater 
than these were set to 225K and 275K respectively. This assumption is valid for the 


images processed in this study since they were from the fall and winter, and it 1s 


unlikely that an emitting surface has a temperature outside this range at that time of 


ear. 


& 


E. SEPARATION ALGORITHM 

Based on the theoretical reflectances discussed above, an automated routine can 
be developed that separates the cloud and surface features. The most important 
separation is between clouds and snow cover and this can be accomplished by 
evaluating the channel 3 reflectance. The lowest theoretical channel 3 reflectance for 
iiquid clouds is 0.11 while for snow the values range from 0.005 to 0.025. There is 
some overlap between ice clouds and snow but a significant number of ice cloud pixels 
can be correctly analyzed along with all liquid clouds. Land can be easily separated 
since its channel | reflectance is about 0.10, well below that of the remaining surfaces. 
When attempting to make a distinction between the remaining ice clouds and snow the 
routine will have to use the channel 3-channel 4 temperature factor. In this case the 
reflective properties are nearly identical in both channels 1 and 3 and it is hoped that 
their different transmissive properties between channels 3 and 4 will improve the 
analysis. 

The ability of an automated routine to separate these surfaces depends on the 
amount of each of them present in a single pixel and this must be considered when the 
separation thresholds are derived. For example, liquid clouds and land covering a pixel 
together may combine to give a reflectance value that is recognized as either liquid 
cloud or land, or possibly neither of these. In the case of clear skies over snow cover, 
the measured reflectance may be significantly less than 0.75 to 0.80 because even with 
deep snow, buildings, roads and trees protrude above the snow cover. This 1s 
especially important in thin snow cover where the underlying surface and adjacent 
areas of no snow affect the measured reflectance. For a single pixel, the measured 


Fenlect ance is expressed py. 


ae 0.709. + 0.60p, + 0.109, + 0.80p, (2.13) 
and 
a 0.259, + 0.05p, + 0.05p; + 0.02p, (2.14) 
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The following examples illustrate how the pixel coverage must be taken into 
account when deciding on the separation thresholds: 


1. If. for example, a decision were to be made between land and thin cirrus based 
Owe cidomel |renectance, Eq. 2:13 becomes: 


r, ~ 0.60p, + 0.109, 


If the threshold channel | reflectance were set at 25 percent, a pixel would have 
to be covered oy at least 30 percent of cloud in order to be recognized as such. 
If the pixel coverage were less, it would be analyzed as land. 


tv 


If a decision were to be made between snow and liquid clouds based on the 
channel 3 reflectance. Eq. 2.14 becomes: 


r; ~ 0.25p, + 0.02p, 


If the threshold channel 3 reflectance were set at 10 percent, the pixel would 
have to be covered by at least 35 percent of cloud in order to be recognized as 
such. If the coverage were less, it would be analyzed as snow. 


Thus when a pixel is covered by more than one of the surfaces to be analvzed, the 
classification that is made depends not only on the surface that dominates but also on 
the combined reflectance of the scene as sensed by the satellite. This is an tmportant 
consideration when one 1s interested in analyzing thin clouds correctly and also in cases 
covered by this study where the difference between snow-covered and snow-free land 1s 
to be analyzed correctly. It is particularly important when using coarse-resolution 
geostationary data where a variable pixel composition is more likely due to the larger 


area covered by the individual pixels. 


Il. DEVELOPMENT OF THE ANALYSIS ALGORITHM 


A. OVERVIEW 

This chapter discusses development of the analysis algorithm. Although the 
primary objective was the separation of low clouds and snow cover, the thresholds 
were developed so the routine would also work well when snow is not present, and in 
the presence of high-level clouds including cirrus. First, observed reflectances from 
channels | and 3 were compared with the theoretical values from Chapter I]. Then the 
observed values of reflectance were used to determine the separation thresholds. Using 
an image analysis work station, land pixels were assigned grayshade values between 0 
and 25, and snow pixels were assigned grayshade values between 26 and 50. This 
darkened the snow pixels in relation to the cloud pixels, and this range was enough to 
show the boundary between snow-free land and snow-covered land, and to show some 
texture within each classification. Clouds were displayed as in an infrared image on a 
grayvshade scale between 75 and 255. The temperature scale 1s 225K to 275K and this 
range allows distinction between low-level warm clouds and high-level cold clouds. 
Clouds were the most important part of the analysis so they were given the greatest 
range of gravshades. However, it is hoped that the land and snow grayshade ranges 
are large enough so various enhancement techniques can be employed to analyze them 


in greater detail if necessary. 


B. SATELLITE DATA DESCRIPTION 

Data for the benchmark cases were taken from the 2000-2100 GMT pass of 
NOAA-I9 on November 9 and 10, 1986 (see Figs. 3.1 through 3.4). A winter storm 
passed through the upper midwestern United States and south-central Canada 
producing the first major snowfall of the 1986-87 season, and both of these passes 
cover that area. The storm left 3-10 inches of snow on the ground in North Dakota 
and South Dakota, and 1-3 inches in Minnesota and western Wisconsin. A typical 
frontal-band cloud pattern is present in the 10 November pass, with low- to mid-level 
enhanced cumulus clouds streaming in behind the frontal system, along with scattered 
snow showers (Fig. 3.3). The images have a combination of clouds over both clear land 
and over snow, and clear conditions with clear land and snow cover. The passes cover 
the area from 40° north latitude to 60° north latitude, and 80° west longitude to 115° 


west longitude. 


30 


Two subscenes were chosen from each pass. each of which covers an area of 
about 500 square kilometers. Subscene | of November 9 (Figs. 3.5 and 3.6) shows 
most of Wisconsin and Lake Michigan, western Michigan and the extreme northern 
portions of Illinois and Indiana. This area was covered by mainly clear land and low 
clouds and was sampled for these surfaces. Subscene 2 of November 9 (Figs. 3.7 and 
3.5) shows the eastern half of Nebraska, western Iowa, and southeastern South 
Dakota. This area was chosen for the cirrus clouds that are present; there was an 
abundance of both thick and thin cirrus here. There was no snow on the ground in 
either of these November 9 subscenes. Subscene | of November 10 (Figs. 3.9 and 3.10) 
has Lake Huron in the center with most of Michigan and the southernmost portion of 
Canada. This area had thick and thin cirrus along with low level clouds available for 
sampling, but no snow on the ground. Subscene 2 of November 10 (Figs. 3.11 and 
3.12) shows most of Minnesota and the eastern portions of North and South Dakota. 
This area had clear skies and snow cover in eastern North and South Dakota and 
western Minnesota, with low-level convective clouds in eastern Minnesota. There were 
a few of these cloud elements over the snow cover in southwestern Minnesota. 

Figs. 3.1 through 3.4 show the November 9 and November 10 passes in what will 
be referred to as the overview format. These figures show channel | of both passes 
along with an image of derived channel 3 reflectance. There are 2048 pixels in each of 
approximately 2500 lines in these passes, but only every eighth line and every eighth 
pixel were processed and displayed. The overviews were then enlarged to cover a full 
512 Sbyesl2 screen.* Data were taken from seven groups of 10 by 10 pixels on the 
November 9 pass and six groups of 10 bv 10 pixels on the November 10 pass, and 
these boxes are labeled in Figs. 3.2 and 3.4. Due to the way the overview for each pass 
was processed, these groups are actually every eighth line and pixel of an 80 by 80 box. 
Care was taken to insure that all of the surfaces to be analyzed were sampled in nearly 
equal numbers, and also that the areas chosen were comprised only of either snow, 
snow-free land, ice clouds or liquid clouds. Uniformity of these areas was determined 
subjectively by evaluating the visible and infrared images along with plotted surface 
observations for the area. From both passes a total of 300 land pixels, 400 ice cloud 
pixels, 300 liquid cloud pixels, and 300 snow pixels were sampled and the following 
information calculated: 

imeeChannel | ¥etlectance. 
2a Channel 3 reflectance. 


3. Channel 4 brightness temperature. 
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Fig. 3.1 9 Nov 86 overview - NOAA-9 AVHRR channel I. 
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Fig. 


Same as Fig. 3.1 except for derived channel 3 reflectance with data samples. 





Fig. 3.4 Same as Fig. 3.3 except for derived channel 3 reflectance with data samples. 
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4. Channel 3-channel 4 temperature factor. 

5. Channel 3 reflectance not including the anisotropic reflectance factor. 
Additionally, the averages and standard deviations of all of these were calculated for 
each 10 by 10 box. 

Figs. 3.5 through 3.12 show channel | of the subscenes along with images of the 
derived channel 3 reflectance. The 10 by 10 boxes of data samples are labeled on the 
images of derived channel 3 reflectance for easy reference to the tables of reported 
values in the following section. For all four subscenes, a total of 600 land pixels, 600 
thin/ warm ice cloud pixels, 600 cold/thick ice cloud pixels, 1000 liquid cloud pixels and 
700 snow pixels were sampled. Separation of the ice clouds in this manner was due to 
the need to more accurately specify thresholds that will separate them from everything 
else. The information calculated was the same as with the overviews except the 
channel 5 reflectance not including the anisotropic factor. This was calculated only for 
the overviews and its significance will be discussed later. The averages and standard 


deviations were calculated for each individual box as was done with the overviews. 


C. DATA SAMPLES 

Tables 3 through 8 show the observed reflectance and the channel 3-channel 4 
temperature factor for the sampled areas from the overviews and subscenes. The range 
of values reported for each box includes one standard deviation either side of the 
average. The solar zenith angle reported with each group 1s for the center of the box. 
Figs. 3.13 and 3.14 summarize these values in graphical form to illustrate how the 
separation thresholds are selected. Theoretical values reported in Chapter II for this 
range of solar zenith angles (64° to 79°) are 0.18-0.42 for liquid clouds, 0.01-0.30 for ice 
clouds, 0.005-0.025 for snow cover and about 0.10 for snow-free land. The major 
differences between theory and observation are with liquid clouds and snow cover. 
With liquid clouds the discrepancy is possibly due to a difference between parameters 
used in model calculations and what was actually observed. The higher values of the 
liquid cloud observations verified quite well. With snow cover the discrepancy is likely 
due to error in the method of calculating the channel 3 reflectance. The assumptions 
affect only the thermal emission portion of the channel 3 radiance measurement which, 
in the case of snow cover, has a magnitude comparable to that of solar reflection. 
Thus a larger error should be expected in the reflectance calculation. These observed 


values verifv that there are sufficient differences on which to base an automated 





Fig. 3.5 9 Nov 86 subscene | - NOAA-9 AVHRR channel I, full resolution image. 
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Fig. 3.6 Same as Fig. 3.5 except for derived channel 3 reflectance with data samples. 
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Fig. 3.7 9 Nov 86 subscene 2 - NOAA-9 AVHRR channel 1, full resolution image. 





Fig. 3.8 Same as Fig. 3.7 except for derived channel 3 reflectance with data samples. 
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analysis routine. The importance of channel 3 in developing this routine is evident 
from Figs. 3.135 and 3.14. Without channel 3, a routine would have to rely on 
distinctions in channel | reflectance and/or the channel 4 temperature. The only 
distinct surface in channel | is land, and the channel 4 temperature could be used onlv 
if the ground temperature were known from a conventional observation. Channel 3 is 
the key to providing the distinction necessary for the automated routine to separate 
these surfaces successfully. 

Fig 3.15 shows the observed channel 3-channel 4 temperature factor for the 
sampled snow cover and cirrus, and there is a great deal of overlap at values greater 
than 20.0. It is likely that most of the cirrus with a lower factor will be analyzed 
correctly because the channel 3 reflectance is high, thus leaving very few to be 
unambiguously separated from snow cover based only on this factor. 

Taylor and Stowe (1984) report the anisotropic reflectance factor for clear land, 
clear snow, ocean and low, middle and high clouds as derived from the NIMBUS-7 
Earth Radiation Budget instrument. These factors are for spectrally-integrated albedo 
between 0.2}tm and 50m. They specify the factor as a function of solar zenith angle 
and what thev refer to as the angular bin, which is determined by the satellite zenith 
angle and the relative azimuth angle. Each bin spans approximately 12-15 degrees of 
Satellite zenith angle and 10-20 degrees of relative azimuth angle, while the solar zenith 
angles are stratified into ten categories from 0 to 90 degrees with an interval of 6 to 25 
degrees. 

Two major assumptions concerning this factor which were applied to this study 
ale: 


1. The factor for the low cloud category adequately describes the factor for all 
clouds. 


to 


Since the factor is used to determine the scene classification, its variation 
among the different surfaces 1s neglected. 


These assumptions were made primarily to ease the burden of computer storage and 
processing. The values of this factor reported by Taylor and Stowe were determined 
with prior knowledge of the scene classification; however, this routine needed the factor 
first in order to make its classification decision. For a particular bin and solar zenith 
angle, the factors for land, snow and clouds are close enough so that the average 1s a 
reasonable estimate for any of the three features. For this reason it was decided to use 
the average anisotropic factor to classify the scene. A single array of factors was used 


in the computer program, with 49 bins and 7 categories of solar zenith angle. The 
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Fig. 3.9 10 Nov 86 subscene I - NOAA-9 AVHRR channel 1, full resolution image. 





Fig. 3.10 Same as Fig. 3.9 except for derived channel 3 reflectance with data samples. 
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Fig. 3.11 10 Nov 86 subscene 2 - NOAA-9 AVHRR channel 1, full resolution image. 
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Fig. 3.12 Same as Fig. 3.11 except for derived channel 3 reflectance with data samples. 


solar zenith angle categories were reduced from 10 to 7 covering the interval from 45° 
to 90°. since all of the images analvzed for this study have the solar zenith angle within 
this interval. With all of these assumptions taken into account, the program chose the 


anisotropic factor from one table of 343 values rather than six tables of 490 values. 


TARDE es 


Observed reflectance for channels 1 and 3 and the channel 3-channel 4 
temperature factor (Nov 9). 


Surface (Group) seth Ch 1 Ch 3 emp Factor Feuthe 
Ice Cloud (1) 76 0.58 - 0.74 0.04-0.06 10.29-15.79  -19 
Ice Cloud”) 79 0.57 - 0:66 O1032 0055017. 63 27 Gt -25 
Ice Cloud (3) 6S 0.64-0.87 0.09-0.27 3.78 - 8.88 15 
Land (4) 67 0.09-0.11 0.03 - 0.04 52.21-113.07  -24 
Land (5) 71 0.13-0.16 0.04-0.06 3836-7400  -29 
Land (6) 71 0.12-0.14 0.05-0.07 40.48 - 60.49 -18 


Liquid Cloud (7) 79 0.63 -0.77 0.14-0.20 9.04 - 11.42 0 


Accounting for anisotropy is important in the calculations and the degree of 
importance can be determined by the amount that the factor is different from 1.0. For 
example, a factor of 0.85 means that if the factor were left out of the calculations the 
reflectance would be about 15 percent less than what it should be. In general, the 
anisotropic factor becomes large when the solar zenith angle, the satellite zenith angle 
and the relative azimuth angle are all large (see Fig. 2.1), NOAA’s polar orbiters have 
both morning and afternoon passes so in all cases the solar zenith angle is large. The 
satellite zenith and relative azimuth angles are large when the sensor is looking in the 
direction of the sun, and this occurred with the morning pass when the sensor looked 
to the east, and for the afternoon pass when the sensor looked to the west. For both 
passes, when the sensor looked towards the sun the relative azimuth decreased 
northward, and when the satellite looked away from the sun, the relative azimuth 


increased northward. There were some instances where the anisotropic factor was 
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Surface (Group) 


Liauid Cloud (1) 
Snow (2) 

Snow (3) 

Snow (4) 

Ice Cloud (5) 
Liquid Cloud (6) 


i 


Liquid Cloud (1) 
Liquid Cloud (2) 
Liquid Cloud (3) 
Liquid Cloud (4) 
Land (5) 

Land (6) 

Land (7) 

Land (8) 

Land (9) 

Land (10) 


: 
Surface (Group) 
| 
| 





Solar 
Zemth 


74 


Solar 
Zemin 
"5 
re 


TABLE 4 


Observed reflectance for channels | and 3 and the channel 3-channel 4 
temperature factor (Nov 10). 


Ch | 


0.39 - 0.60 
0.49 - 0.61 
Oe 0055 
0.53 - 0.69 
0.60 - 0.73 
0.48 - 0.67 


Ch 3 


0.20 - 0.36 
0.02 - 0.03 
0.03 - 0.04 
0.02 - 0.03 
0.02 - 0.04 
Re ae 0 ae 


DSU ee sae) 


Observed reflectance for channels 1 and 3 and the channel 3-channel 4 
temperature factor (Nov 9-subscene 1). 


Ciel 


0.43 - 0.55 
0.41 - 0.53 
0.64 - 0.72 
0.60 - 0.67 
0.14 - 0.17 
0.14 -Q.17 
0.16 - 0.18 
Odio 0.15 
Olo= 0.17 
Oro 0.18 


4} 


1.77 - 10.83 
20.09 - 23.7] 
20 = 35.06 
19:61 - 29.87 
10.53 - 14.98 
Bog = 20.9) 


temo Pacser Pct Che 


Ch 3 


0.16 - 0.19 
0.16 - 0.20 
Oe a0 
0.13 - 0.16 
0.04 - 0.06 
0.04 - 0.06 
0.05 - 0.07 
0.04 - 0.05 
0.05 - 0.07 
KE05 - 0710 


Temp factor 


8.58 - 9.50 

8.46 - 9.63 

8.01 - 8.97 

HUGE WU OL| 
59° 33'=365,05 
43.21 - 61.93 
34.31 - 46.74 
ee it 00,9) 
Spey pues! os Wi 
35,07 35.19 





TABI SG 


Observed reflectance for channels | and 3 and the channel 3-channel 4 
temperature factor (Nov 9-subscene 2). 


Surface (Group) Pea lial Cho. Temp Factor 

Ice Cloud (1) 65 0:47 Saino 0.06 - 0.08 523 26,05 

| Ice Cloud (2) 65 Oe Oa 0.06 - 0.08 ee SKE 

Ice Cloud (3) 64 0.57 - 0.63 0.06 - 0.07 5.29 - 6.05 
Ice Cloud (4) 67 0.26 - 0.29 0.08 - 0.09 12.38 - 14.02 

Ice Cloud (5) 68 0.19 Gri 0.06 - 0.07 21.92 - 24.64 
Ice Cloud (6) 67 O15 00kd 0.04 - 0.06 30.59 - 39.37 





TAB Ea 


temperature factor (Nov 10-subscene 1). 


Temp Factor 


Ice Cloud (1) 
Ice Cloud (2) 
Ice Clouda, 
Ice Cloud (4) 
Kee @loudss: 


Liquid Cloud (6) 
Liquid Cloud (7) 
Liquid Cloud (8) 


Surface (Group) aan Ch 1 meen 3 
76 0.34 - 0.44 0.03 - 0.05 
76 0.34 - 0.45 0.02 - 0.04 
74 0.53 - 0.56 0.02 - 0.03 
US 0.61 - 0.68 0.05 - 0.06 
75 0.53 - 0.62 0.05 - 0.06 
76 0.55 - 0.58 0.08 - 0.09 
76 0.50 - 0.62 0.09 - 0.13 
76 0.47 - 0.62 0.10 - 0.16 
fo 0.44 - 0.61 0.09 - 0.18 


ice Cloud) 


Observed reflectance for channels | and 3 and the channel! 3-channel 4 


25.86 - 34.82 
A D.22 
De oii O es UB as) 
9 OOS EL 
9.06 - 11.00 
18.22 - 19.46 
13.75 - 17.45 
12.33 - 15.84 
5.56 - 7.44 
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Observed reflectance for channels | and 3 and the channel 3-channel 4 
temperature factor (Nov 10-subscene 2). 


Solar 
Surface (Group) Zenith 


Chal Ch) Temp Factor 
Snow (1) 7] 0.68 - 0.76 0.03-0.04 22.98 - 28.53 
Snow (2) 7 0.77 - 0.82 0.03-0.04 23.87 - 33.15 
Snow (3) 7 0.75 - 0.81 0.03-0.04 24,92 - 32.00 
Snow (4) 70 0.82 - 0.83 0.02-0.03 23.14 - 28.89 
Snow (5) 7] 0.81 - 0.85 0.02-0.03 26.84 - 34.64 
Snow (6) 70 0.76 - 0.81 0.03-0.04 23.14 - 28.89 
Snow (7) 70 0.57 - 0.69 0.02-0.03 32.57 - 41.51 
Liquid Cloud (8) 72 0.42 - 0.56 0.21 - 0.30 4.51-8.18 | 
Liquid Cloud (9) 7! 0.43 - 0.62 0.23 - 0.34 208 = Un 
Liquid Cloud (10) 71 0.42 - 0.63 0.19 - 0.33 y7eote 





significantly less than 1.0. This occurred at various combinations of sun-satellite 
geometry when the satellite sensor was not looking towards the sun. 

Tables 3 and 4 show the change in the channel 3 calculated reflectance, as a 
percentage of the channel 3 reflectance, neglecting the anisotropic reflectance factor. A 
negative change means that the reflectance is less when the factor is neglected. 
Referring to Eq. 2.12, a negative change occurs when the factor is less than 1.0. so the 
denominator of Eq. 2.12 is larger when the factor is neglected. As discussed in 
Chapter II. when the factor is less than 1.0 the calculated reflectance would vield an 
underestimate of the radiant exitance if anisotropy were neglected, and the calculated 
reflectance would yield an overestimate when the factor is greater than 1.0. 

A closer analysis reveals that the anisotropic factors ranged from a minimum of 
0.59 to a maximum of 3.41, meaning that there was as much as a four-fold variability 
in reflectance due to anisotropy alone, clearly indicating that this eflect cannot be 
neglected. The factors for these surfaces in the tables of Taylor and Stowe show values 


as high as 5.98 but the values applicable to this study were restricted to this range 
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Fig. 3.13. Summary of NOAA-9 AVHRR observed channel | reflectance 
for 900 Jand pixels, 1600 ice cloud pixels, 
1300 liquid cloud pixels, and 1000 snow pixels. 


because the satellite zenith angle at the edge of these passes was approximately 69°. 
Thus the angular bins were restricted to between I and 4]. This fact can be used to 
effect a slight improvement in computer processing. Fewer data values to be stored 
means the program will run slightly faster and that is important with this 
computationally-intensive routine that chooses a factor for a total of 262.144 pixels in 
a single image. There are some instances where the factors for the various scenes varv 
Widely within a particular bin and solar zenith angie category. In these cases the 
average 1s not a very good estimate for the program to use and it may be better to 
sacrifice an improvement in speed to gain more in accuracy of the anisotropic factor. 
Now that fewer values need to be stored, instead of using an average for all scenes, a 
table of values for each individual tvpe of classification could be used in the program. 
A preliminary guess of the scene classification could be made and then an iterative 


process used to arrive at the correct factor for that particular classification. 
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Fig. 3.14 Same as Fig. 3.13 except for derived channel 3 reflectance. 


The following are important points from both overviews concerning the 
anisotropic factor: 


1. November 9 (Table 3): For Group 7 the satellite was looking away from the 
sun so negiecting amisotropy had a verv small effect on the calculated 
reflectance. Contrast this with Group 3 where the satellite was looking towards 
the sun. In this case neglecting anisotropy had a significant effect. For Group 
7 tne relative azimuth was about 25-26° and the amisotropic factor was constant 
at 1.0. For Group 3. the relative azimuth was about 134-136" and the 
anisotropic factor varied between 0.96 and 1.22. making a positive change in 
reiiccmince @ehen amsotropm was neglected. It 1s mteresting to note that the 
change for land was verv significant. Groups 4, 5 and 6 are relatively close to 
the satellite subpoint but the anisotropic factors all came out to be between 
0.75 and 0.85 and this translated into a large percentage change. This may be 
due to the fact that there was a large difference between the average factor used 
in the program and the factors for land reported by Tavlor and Stowe. 


2. November 10 (Table 4): In this image, both Groups | and 6 had a small 
change. For Group 6 the relative azimuth was very small. For Group | the 
relative azimuth was larger than for Group 6, but the solar zenith and satellite 
zenith did not work together to make the factor very much different from 1.0. 
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Fig. 3.15 Same as Fig. 3.13 except for channel 3-channel 4 temperature factor, 
1000 snow pixels and 1600 ice cloud pixels. 


Group 5 is near the satellite subpoint and the factor here was constant at 0.74, 
giving a large percentage change when amisotropy is neglected. The most 


a 


interesting values for this image are for the snow cover of Groups 2, 3 and 4. 
These groups are fairlv close together and the sun-satellite geometry was similar 
for them all, however it was different enough to cause a quite variable change in 
reflectance among the 3 groups when anisotropy 1s neglected. 


To summarize, accounting for anisotropy is important when calculating 
reflectance. The change in reflectance becomes important (greater than 10 percent) 
when the factor 1s greater than approximately 1.10 or less than approximately 0.90. It 
Seems as though the greatest effect of neglecting anisotropy should be when the 
satellite is looking towards the sun, but it was shown here that for these images it 1s 
also important near the Satellite subpoint with a large solar zenith angle. With the 
former the effect is large due to forward scattering, while with the latter the effect 1s 
large due to the small degree of backward scattering. The factors used here are 


admittedly only first-order approximations and in some instances provide only verv 
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crude estimates. Taylor and Stowe state that some of the angular bins were very 
sparselv populated and subject to large error in the estimate for that bin. The factors 
could be improved by making the angular resolution of the bins finer and taking more 
observations for the sparsely populated bins. Measurements should also be taken in 


smaller spectral intervals to address the possibility of spectral variability of the factors. 


D. ALGORITHM DEVELOPMENT 

The algorithm was developed in phases, and additional data were gathered as 
necessary to determine the conditions under which the algorithm failed. Analyses for 
each phase were evaluated subjectively and compared with visual and infrared imagery 
along with plotted surface observations for the area. Subscene 2 of November 10 1s 
the only one of these images with snow cover. Since the focus of this studv is the 
problem of snow,low cloud discrimination, this image will be included again in Chapter 
TV where it will be evaluated against the individual surface observations in the area and 
its verification statistics included with the other images evaluated there. 

In order to satisfy the primary objective of this analysis routine, 1t was decided 
that there should be three steps executed in order during the classification process: 


1. Separate liquid clouds and as many of the ice clouds as possible from 
Every ining else. 


2. Separate land from snow and the remaining ice clouds. 

3. Finally, separate snow and the remaining ice clouds. 
This approach was adopted because the distinction between clouds and everything else 
is the sharpest. Anv ambiguities that arose between low clouds and snow or land were 
resolved as clouds since, for operational forecasting, it 1s better to overestimate the low 
cloud cover than to underestimate it. This logic was employed to correctly separate all 
snow and low clouds, then anv confusion between snow and ice clouds could be 
resolved later by manual analysis of the infrared imagery. 

All images of the analyses presented in this section emphasize the cloud cover 
analvsis. First the pixel classification was made, then land pixels were assigned a 
gravshade value between 0 and 25 and snow pixels were assigned a grayshade value 
between 26 and 50. The actual value was assigned based on the channel | reflectance. 
All snow and snow-free land areas appear dark in these images but enough variation 1s 
actually present to analyze the surface features in greater detail by using an image 
enhancement technique. The cloud pixels were assigned a grayshade value from 75 to 
255 based on their channel 4 temperature. This provides some texture to the cloud 


analvsis and also shows variability in cloud height. 
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According to Fig. 3.14, liquid clouds and most ice clouds can be separated 
from land and snow by specifying a threshold with the channel 3 reflectance alone. A 
reasonable choice for the boundary between liquid clouds and snow/land is somewhere 
between 7 and 10 percent. Snow should be separated well but there is some overlap 
between clear land and clouds since the channel 3 reflectance of land was as high as 10 
percent. There are verv few land pixels with a reflectance this high, so a threshold of 
7.7 percent was decided upon for this first attempt. Ice cloud reflectance ranged from 
2 to 27 percent so many of the ice clouds will make this separation. 

Channel 3 reflectance cannot be used to separate land from snow cover and 
the remaining ice clouds (see Fig. 3.14). There is adequate distinction in the channel | 
reflectance (Fig. 3.13), and a threshold of 20 percent was decided upon to accomplish 
this separation. This allows thin snow cover and thin ice clouds to be analvzed 
correctly if the channel | reflectance 1s at least 20 percent, thus making the pixel 
approximately 25-30 percent cloud-filled. If snow or ice clouds have a channel 1 
reflectance lower than 20 percent and are thus incorrectly analyzed as land, then the 
clouds are probably so thin that it will not greatly affect the accuracy of the analysis. 

The final separation is between snow and the remaining ice clouds. Neither 
the channel 1 nor channel 3 reflectance can be used to distinguish them, so this 1s 
where the channel 3-channel 4 temperature factor was used (Eq. 2.7). This factor was 
generally greater than 20 for snow cover, although there was some overlap with ice 
clouds (Fig. 3.15). For phase 1 the threshold for this factor was set at 20. 

Summarizing the logic for the analysis routine at this phase of development: 


1. Pixels with a channel 3 reflectancesencatemtmansor caualaton).) eercemmvcre 
analyzed as cloud. Land and snow/ice clouds were passed on. 


tv 


Pixels with a channel 1 reflectance less than 20 percent were analyzed as land. 
Snow/ice cloud pixels were passed on. 
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Pixeis with a channel 3-channel 4 temperature factor of greater than or equal to 
20 were analyzed as snow and the remaining pixels were analyzed as cloud. 
2. Results of Phase 1 Test 
Figs. 3.16 and 3.17 show the results of applying the algorithm to the 
November 9 subscenes. For subscene I, many of the clear land pixels in Wisconsin 
and northern [Illinois were analyzed as cloud. According to the additional data samples 


that were taken from this area, the misanalysis occurred: because the channel 3 
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reflectance of some land pixels was as high as 10 percent. The liquid clouds were 
handled very well but much of the cirrus north of Lake Michigan was analvzed as 
snow because it did not meet the criteria for the first separation, and then in the final 
separation the temperature factor was between 30 and 40, well within the range for 
snow. For subscene 2, practically all of the cirrus was handled correctly, including 
much of the thin cirrus. Some of this was due to the clouds making the reflectance 
threshold, but much was due to the temperature factor. As noted in Table 6, only a 
few of these pixels had a channel 5 reflectance greater than 7.7 percent, but many had 
a temperature factor less than 20. When the channel I reflectance of the thin cirrus 
went below 20 percent it was analyzed as land, as expected. The remaining cirrus that 
was missed was due to its being greater than the threshold for the temperature factor. 

Figs. 3.18 and 3.19 show the results of applying the algorithm to the 
November 10 subscenes. On subscene 1, the difficulty with cirrus is evidenced by the 
gaps in the cirrus band between Lakes Michigan and Superior. The routine handled 
the low clouds over Lake Huron well, but there are some gaps indicating a channel 3 
reflectance that is lower than 7.7 percent. Land was analyzed correctly in this image. 
For subscene 2, it did an outstanding job in separating the snow and low clouds. Most 
of the convective clouds in this image can be picked out on the visual image (Fig. 3.11) 
due to their texture difference with the smooth snow cover. This analysis maintained 
the texture and structure of these cloud streets, and it is easy to see where the snow 
depth increases. Snow pixel grayshades range from 26 to 50 according to their channel 
1 reflectance, thus the deeper snow (greater reflectance) is displayed on the higher end 
of this range. Additional data were taken from eastern North Dakota to further verifv 
Ene ‘charinel @amemectance threshold of 7.7 percent. Ime highest measured channel 3 
reflectance of the snow in this area was 5.4 percent. 

3. Phase 2 

In the first separation, too many land pixels were mistakenly analyzed as 
clouds because the channel 3 reflectance was higher than 7.7 percent. Also, many of 
the ice clouds were incorrectly analyzed as snow. The problem with the land pixels was 
handled by adding a threshold for the channel | reflectance of 20 percent. At the same 
time there was plenty of room to adjust the channel 3 reflectance threshold downward. 
ft was adjusted to 6.5 percent to allow more ice cloud pixels to be analyzed correctly 
and not interfere with the excellent results in the snow area obtained on subscene 2 of 


November 10. An adjustment was also made in the temperature factor. The 
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Fig. 3.19 Same as Fig. 3.16 for 10 Nov 86 subscene 2. 
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additional data taken from the snow cover of eastern North Dakota showed that the 
temperature factor was as low as 15.7. It was understood that adjusting this threshold 
downward would worsen the analysis of cirrus but 1t was lowered to 15.0 to ensure that 
all the snow cover was analyzed correctly. 

Summarizing the logic for the analysis routine at this phase of development: 


1. Pixels with a channel 3 reflectance of greater than or equal to 6.5 percent AND 
channel 1 reflectance of greater than or equal to 20 percent were analyzed as 
cloud. Land and snow/ice cloud pixels were passed on. 


2. Pixels with a channel | reflectance of less than 20 percent were analyzed as 
land. Snow/ice cloud pixels were passed on. 


3. Pixels with a channel! 3-channe! 4 temperature factor of greater than or equa! to 
15 were analyzed as snow and the remaining pixels were analyzed as cloud. 


4. Results of Phase 2 Test 

Figs. 3.20 and 3.21 show the results of application to the subscenes of 
November 9. These adjustments completely solved the problem of land pixels being 
analvzed as cloud but the analysis of ice clouds was worse than in phase 1. Apparently 
lowering the temperature factor to 15 sacrificed more of the ice cloud pixels than 
lowering the channel 3 reflectance added. On subscene 2, the low clouds in Iowa, 
South Dakota and Minnesota were handled worse than in phase 1. The additional 
samples taken from this area showed that the clouds that were missed met neither the 
channel 1 nor the channel 3 thresholds and were so thin that they were analyzed as 
land. 

Figs. 3.22 and 3.23 show the results of application to the subscenes of 
November 10. On subscene | it seems again as, though the adjustment in the 
temperature factor overwhelmed what should have been an improvement in the ice 
cloud analysis due to lowering the channel 3 reflectance threshold, while on subscene 2 
the excellence of the snow/cloud analysis was maintained. 

5. Phase 3 

In order to improve the low cloud analysis in Iowa, South Dakota and 
Minnesota on subscene 2 of November 9, the channel | reflectance threshold was 
lowered to 19 percent. Overall, more ice clouds should also be analvzed correctly as a 
result of this change. The highest observed reflectance for land was 20.3 percent. This 
value was observed in one pixel on subscene | of November 9, and out of 900 pixels 


from the original sample, only 3 had a channel | reflectance between 19 and 20.3 
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Fig. 3.20 9 Nov 86 subscene 1 - phase 2 cloud cover analysts results. 





Pereeeereneedie as Fie. 3.20 tor 9 Nov 86 Subscene 2. 
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Fig. 3.23. Same as Fig. 3v0sor lO Nowe snbscence 


percent. Lowering the threshold will cause these pixels to be analyzed as cloud since 
their channel 3 reflectance also meets the threshold for cloud, but in this particular 
subscene, the low cloud will be overestimated by much less than 1 percent of the total 
land pixels. The tradeoff was an overall improvement in the cloud analysis by correctly 
analvzing more low clouds and ice clouds. 

The highest measured channel 3 reflectance for snow was 5.4 percent. Thus is 
approximately 1.5 percent higher than theory suggests the maximum ought to be so it 
was assumed that the channel 3 reflectance would be no higher than this. The channel 
3 reflectance was lowered in this phase to 5.7 percent 1n a final effort to get more ice 
clouds analvzed correctly without interfering with the accuracy of the snow cover 
analysis. The temperature factor was not adjusted further. Summanizing the logic for 
the analvsis routine at this phase of development: 


1. Pixels with a channel 3 reflectance of greater than or equal to 5.7 percent AND 
channel | reflectance of greater than or equal to 19 percent were analyzed as 
cloud. Land and snow/ice cloud pixels were passed on. 


tr 


Pixels with a channel 1 reflectance of less than 19 percent were analyzed as 
land. Snow ice cloud pixels were passed on. 


3. Pixels with a channel 3-channel 4 temperature factor of greater than or equal to 
15 were analyzed as snow and the remaining pixels were analyzed as cloud. 


6. Results of Phase 3 Test 

Figs. 3.24 and 3.25 show the results of applying the algorithm to the 
November 9 subscenes. All clouds were handied better as a result of the changes made 
in this phase. As expected, a verv few of the land pixels were analyzed as cloud, 
primarily in Michigan and Indiana, but this did not significantly affect the accuracy of 
the analysis. There was an improvement in the ice cloud analvsis but still there was a 
significant portion incorrectly analyzed as snow. This was evident particularly in 
subscene 1. In this image it 1s possible to claim only about a 40 percent accuracy in 
the analysis of the ice clouds, however for subscene 2 an accuracy of about 70 percent 
can be claimed. With subscene 2, a particularly good analysis was obtained with the 
thick cirrus in Nebraska because with these clouds the temperature factor 1s between 5 
and 10. Even as this cirrus thins out over Iowa the analysis remained good until finally 
the channel | reflectance dropped below 19 percent and it was analyzed as land. This 
deficiency in the analvsis routine can be resolved by manual analysis of the coincident 


visual and infrared imagerv. 


Figs. 3.26 and 3.27 show the results of applying the algorithm to the 
November 10 subscenes. The changes made in this phase did improve the cloud 
analysis overall for subscene 1. About 50 percent of the ice clouds were analyzed 
correctly while nearly all the low clouds were handled correctly. On subscene 2, the 


analysis of the low clouds and snow cover was not affected by these changes. 


E. SUMMARY 

Fig. 3.28 summarizes in graphical form the reflectance thresholds that were 
decided upon from the experiments just discussed. The dotted line at a channel 3 
reflectance of I5 percent for land indicates the maximum expected channel 3 
reflectance although the algorithm actually allows the value to be higher. Ice clouds 
Were separated from snow cover using the channel 3-channel 4 temperature factor 
threshold of less than 15.0. As suspected, this did not add much to the analvsis of 
cirrus because solar reflection added too much variability when trving to make a 
distinction between cirrus and snow cover. About 50 percent of the cirrus clouds were 
handled correctly with the channel 3 reflectance threshold of 5.7 percent, and this 
factor added approximately 2 to 5 percent to that. The variability of reflectance for 
both surfaces below 5.7 percent was enough to make unambiguous separation of them 
possible for only a very small range of optical depths. Very few of the remaining ice 
clouds after the first separation were within this range. 

Assuming the data samples are representative of the surfaces being analyzed, 
enough data were taken from the images so that these thresholds should be applicable 
in general. Reflectance is a conservative property of these surfaces so the only 
significant restriction on use of this routine is the presence of daylight. This effectively 
eliminates polar regions during winter, but where there is sunlight this routine should 
work Well in general. This will be tested in the next chapter, where the routine will be 
applied to another group of images with different solar illumination and under varying 


conditions of cloud cover and land/snow cover. 
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Pic ome odie as Fig. 5-24 for 9 Nov 86 supscene 2. 
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Fig. 3.27 Same as Fig. 3.24 for 10 Nov 86 subscene 2. 
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Fig. 3.28 Reflectance thresholds for the automated analysis routine. 


IV. APPLICATION OF THE ANALYSIS ROUTINE 


A. OVERVIEW 

The performance of the analvsis technique developed for this study is evaluated 
in this chapter. Data from two different NOAA satellites, NOAA-9 and NOAA-10, 
were analyzed to add independence to the statistical results. A total of 6 images are 
presented, all of which are independent of each other in terms of areal coverage or 
different satellite. Three of these images were taken from NOAA-9 and three from 
NOAA-10. Two images were analyzed from the November 10, 1986, pass discussed in 
Chapter III, including subscene 2 of November 10. Four additional images were taken 
from a data set from the period January 17-23, 1987. A series of winter storms had 
recently passed through the central United States leaving significant snow cover from 
Oklahoma and Kansas northeastward through Missouri, lowa, and Illinois. The data 
were processed as described in Chapter III, but onlv the subscenes are presented. Each 
of the analvses were evaluated subjectively and then a statistical evaluation 
accomplished with the available surface observations. Finally, the statistical results 
from these images were totalled to assess the general performance of the analvsis 


routine. 


B. STATISTICAL EVALUATION PLAN 
It was assumed that the ground observer can see approximatelv 30 kilometers, so 

the reporting stations were located on each image’ and each pixel within about 30 
kilometers was sampled for cloud cover. The number of pixels sampled was either 709 
or 1257 depending on the resolution for the particular image. The sampled area was 
drawn around each of the reporting stations in the images shown here. Two of the 
images were taken close to the edge of the pass so a [5-pixel radius was sampled for 
these. The other four were taken closer to the satellite subpoint where the resolution is 
better and for these a 20-pixel radius was sampled. All the cloud pixels from these 
samples were then tallied for each station and the cloud cover was determined bv the 
ratio of the number of cloud pixels to the total number of pixels sampled. The percent 
coverage was then converted to four categories as follows: 

1. Clear - less than 2 percent coverage. 

2, scattered - 2 to 50 percent covenase: 


3. Broken - 50 to 98 percent coverage. 
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4. Overcast - greater than 98 percent coverage. 
The total skv cover from the surface observation was compared with that calculated 
from the satellite analysis. 

The statistics for each image are presented in three different groups. The first 
group of statistics reported for each image is those for which both observations either 
side of the analvsis were the same. Included in this group are those stations where 
only one observation was available, whether or not it was the closest observation to 
the time of the analysis. [Included in the second group of statistics is those stations for 
which there was a one-categorv difference in the two surface observations. The 
analvsis was considered to have matched the observation if there was a match between 
the analvsis and either of the observations. If there was no match then the difference 
was evaluated based on the closest available observation. Finally, the third group 1s 
those stations for which there was a two- or three-category difference in the two 
surface observations. In this group, the analysis was considered to have matched the 
observation if it matched either observation or was any category in between the two 
observations. This happened infrequently but occurred when the sky cover changed 
rapidly during the hour. It can be argued that since the analysis took place in between 
the observations, the analvsis 1s accurate if it 1s any category in between or if it 
matches either observation. The analysis was scored as incorrect if, for example, the 
observations were broken and clear with the analvsis showing overcast for that station. 
As with the second group, if the analvsis was incorrect the difference was based on the 
closest available observation. For all groups, the reported percentages include the 
stauons where the analysis and observation matched exactly, and those stations for 
which there was a one-category, a two-category, and a three-category difference 
between the observation and the analysis. 

This approach was taken because it was not possible to verify the satellite 
analvsis with a surface observation taken at the exact same time. Surface observations 
are available on the hour and there was a time difference of 10 to 50 minutes between 
these observations and the time of the satellite image. Most of the observations for the 
hours either side of the satellite pass were the same, leaving little doubt the same skv 
cover was also present at the analysis ttme. There were many cases, however, where 
the observations either side of the analvsis were different, making it impossible to know 


the exact skv cover at analvsis time. 
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C. RESULTS OF THE ANALYSES 

Two images on each of three different days were analyzed and evaluated. For 
each dav, each image was independent of the other based on either different areal 
coverage or different satellite. These cases are described in this section along with the 
results of the analyses. The images shown in this section are the analyses with the 
location of the reporting stations shown. Refer to Appendix C for comparison of these 
analyses with the coincident visible and infrared images. 

1. November 160, 1986 

These first two cases to be evaluated came from NOAA-9I, 2035 GMT, passes 
of the upper midwestern portion of the United States. The first major storm and 
snowfall of the 1986-87 season had occurred in this section of the country, leaving 6-10 
inches of snow on the ground. The cloud pattern was typical, with a weil defined band 
of jet stream cirrus and a large area of enhanced low-level cumulus streaming into this 
area as the storm progressed eastward. A complete sequence of channels 1, 3 and 4 
images along with the analyses of these cases 1s shown in Appendix C, Figs. C.13-C.16 
(Case 1) and Figs. C.17-C.20 (Case 2). 

Case | shows eastern North and South Dakota, western Minnesota and the 
extreme southern portion of Quebec, Canada. The results of this analysis are shown in 
Fig. 4.1. The deepest snow cover is in northern North Dakota and Minnesota and the 
depth decreases southward. The clouds in this image are all low-level enhanced 
cumulus with bases at about 1500 meters and are approximately 100 to 500 meters 
thick. Table 9 lists the observed sky cover for 2000 GMT and 2100 GMT, along with 
the analyzed sky cover. The results of the statistical analysis are summarized in Table 
10. Case 2 is southwestward of the area shown in Case 1 and shows southeastern 
South Dakota, southwestern Minnesota, eastern Nebraska and western Iowa. The 
results of this analysis are shown in Fig. 4.2. The snow cover in this image is rather 
thin for the most part, but increases in depth in western Nebraska and South Dakota. 
The cloud cover 1s a portion of the frontal band and jet stream cirrus, with a small 
portion of low-level cumulus in central South Dakota and southwestern Minnesota. 
Table 11 lists the observed sky cover for 2000 GMT and 2100 GMT, along with the 
analyzed sky cover. The results of the statistical analysis are summarized in Table 12. 

The routine produced excellent results with these cases. For Case | there is a 
definite separation between the low clouds on the right hand side of the image and the 


snow cover on the left hand side, and there is some low cloud over snow in the center 
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Fig. 4.1 Case 1 cloud cover analvsis results. The location of the 
surface reporting stations is at the center of the labeled circles. 


(see Appendix C, Figs. C.13-C.15). The general pattern evident from the visible and 
infrared images 1s apparent in the analysis. In both cases, the section of clouds over 
snow cover in Minnesota and South Dakota can be seen from the texture differences in 
the visible image and this, along with the entire snow covered area. was handled very 
Well [he cirrus shield in Case 2 was handled generally well but, as expected, the 
analvsis had difficulty with the thinner cirrus in central Jowa. 

At FAR (Fargo, ND) on Case 1, there are some snow pixels that had a 
channel 3 reflectance above the 0.057 threshold, but there also appear to be some 
genuine cloud pixels within the area sampled. However these cloud pixels are so few 
Maretneye were POSsIOly mot seen by the observer. In Case 2) three stations were 
affected by the routine’s difficulty with cirrus; DSM (Des Moines, IA). FRM 
(Fairmont, MN), and OTG (Worthington, MN). In this area, the channel 3 
reflectance of the cirrus was below threshold and the channel 3-channel 4 temperature 
factor did not help appreciably. The factor did, however, add to the cirrus analysis in 


southwestern lLowa. 
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TABLES 


Analyzed versus reported sky cover for Case ]. 


Reported Analyzed Reported 
Station Coverage-20Z Coverage Coverage-21Z 
olen Clear Clear Clear 
AXN Overcast Scattered Overcast 
BJ Scattered Scattered Scattered 
BRD Scattered Scattered Not available 
DVL Overcast Scattered Not available 
FAR Clear Scattered Not available 
GFK Clear Clear Gear 
Hie Scattered Broken Scattered 
INL Scattered Broken Not available 
IMS Clear Clear Clear 
Se Scattered Broken Scattered 
Te Clear Clear Overcast 
Yon Scattered Broken Scattered 


There are some discrepancies between the analysis and the observations due 
perhaps to differences in perspectives of the satellite and observer. Two particularly 
close calis in Case 1 are STC (St. Cloud, MN) and INL (International Falls, MN). 
The sky cover for both stations was composed of individual cumulus elements scattered 
around the sky and it 1s difficult for the observer to evaluate the total sky cover under 
these conditions. The analyzed skv cover was 70 percent for INL and 62 percent for 
STC which 1s close enough to the observed scattered condition to make this a possible 
explanation. In Case 2, PIR (Pierre, SD) and RWF (Redwood Falls, MN) had an 
analvzed sky cover that was much less than observed. This also occurred at DVL 
(Devils Lake, ND) on Case 1. At all three of these stations the cloud cover appears to 
be directly over the airfield, perhaps leading the observer to evaluate it as overcast even 
though the horizontal extent of the cloud was not great. The remaining discrepancies 
for Case | were at AXN (Alexandria, MN), HIB (Hibbing, MN), and YQK (Kenora, 
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Verification of Case 1 cloud cover analysis results. 


Category Count Percentage 
Group | Correct ay le 42 
l-category difference SO 42 
2-category difference Ze 16 
3-categorv difference Ourz 0 
Group 2 Correct 0/0 0 
l-category difference 0/0 0 
2-category difference 0/0 0 
3-category difference 0/0 0 
Group 3 Correct 1/1 100 
l-category difference OI 0 
2-category difference Ol 0 
3-category difference 0/1 
Totals Correct GAs 46 
l-categorv difference Sis 38 
2-categorv difference 25 16 
3-category difference 0/13 0 





Quebec, Canada). It is difficult to attribute these to a weak spot in the analysis 
algorithm; they may be due to a difference in perspectives of the observer and the 
satellite. For Case 2 the remaining discrepancies were at FSD (Sioux Falls, SD), OFK 
(Norfolk, NE), and OLUL (Columbus, NE). At both OLU and OFK the problem with 
cirrus again was the culprit and the temperature factor did not add very much to the 
cirrus analysis. At OLU the analyzed sky cover was 95 percent, which is within 3 
percent of the observed sky cover. At FSD, all analyzed brightness temperatures were 
within 3 degrees of the observed shelter temperature, so any of the scattered clouds 
that were reported are either below the resolution of the satellite or too thin to be 


analyzed correctly. 
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Fig. 4.2 Same as Fig. 4.1 for Case 2. 


2. January 17, 1987 

The next two cases to be evaluated were from the same area but from different 
satellites and were separated by about 6-7 hours. Major storms the previous week 
deposited a very deep layer of snow and the snow depth is about 10-16 inches in 
eastern Colorado, thinning out in western Kansas. There is a frontal cloud band and 
cirrus shield in these images and some low-level clouds in the cold air to the north. A 
complete sequence of channels 1, 3, and 4 images along with analvses for these cases 1s 
shown in Appendix C, Figs. C.21-C.24 (Case 3) and Figs. C.25-C.28 (Case 4). 

Case 3 shows eastern Colorado, southwestern Nebraska, western Kansas and 
the Oklahoma panhandle. The deep snow in eastern Colorado thins out and seems to 
merge with an area of low clouds in southwestern Nebraska. The frontal band and 
cirrus shield extend northeastward from northern New Mexico into south-central 
Nebraska. This image was taken from the NOAA-9 pass of 2135 GMT and the results 
of the analysis are shown in Fig. 4.3. Table 13 lists the observed sky cover for 2100 
GMT and 2200 GMT, along with the analyzed sky cover, and the results of the 
analysis are summarized in Table 14. Case 4 shows generally the same area as Case 3 


but a little farther to the south. This image shows more of the frontal band and cirrus 
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Analyzed versus reported sky cover for Case 2. 


Reported Analyzed Reported 
Station Coverage-20Z Coverage Coverage-21Z 
ose Overcast Overcast Overcast 
ANW Clear Clear Clear 
BBW Clear Clear Clear 
BRX Not Available Clear Clear 
DSM Broken Scattered Not available 
EAR Overcast Clear Clear 
FOD Overcast Overcast Broken 
FRM Overcast Broken Not available 
feo 1) Scattered Clear Not available 
GRI Broken Broken Overcast 
HON Clear Clear Clear 
EBE Clear Clear Clear 
NK Broken Broken Broken 
Niel Clear Clear Clear 
MHS Clear Clear Clear 
OF}. Overcast Broken Not available 

OLUL Overcast Broken Overcast 
OMA Overcast Overcast Overcast 
OG Not available Broken Overcast 
PIR Overcast Scattered Overcast 
RWF Overcast Scattered Not available 

| Sex Overcast Overcast Broken 


VIN 


Clear 


Clear 


Clear 
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Verification of Case 2 cloud cover analysis results. 


Category Count Percentage 
Correct 11/19 58 
l-category difference 6/18 
2-category difference 2/19 


3-category difference 0/19 


Connect 3/3 
l-category difference 0/3 
2-category difference 0/3 
3-categorv difference 0/3 

Correct 1/1 
l-category difference 0/1 
2-category difference 0/1 
3-categorv difference 0/1 

Totals Correct 15/23 
l-category difference 6/23 
2-category difference 2128 


3-category difference 0/23 





shield and there are low clouds present over the snow cover in eastern Colorado and 
western Kansas. Additionally, there is a widespread area of low- to mid-level clouds in 
the Texas panhandle. The results of the analysis are shown in Fig. 4.4. Table 15 lists 
the observed sky cover for 1500 GMT and 1600 GMT, along with the analyzed skv 
cover. This image came from the NOAA-10 pass of 1513 GMT and the results are 
summarized in Table 16. 

In spite of the fact that there is a great deal of cirrus in these images, the 
analysis routine worked quite well. Particularly encouraging is the analvsis of the low 
clouds and snow, which was virtually all correct. The most interesting point about 
these analyses is the role of shadows cast due to the large solar zenith angle. This was 


not considered at the beginning as a possible source of error but these analyses 
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Fig. 4.3 Same as Fig. 4.1 for Case 3. 


illustrate the effect quite graphically. Case 3 was an afternoon pass and the sun was to 
the southwest casting shadows on the snow cover in Colorado and Kansas. Case 4 
was a morning pass and the sun was to the southeast casting shadows on the snow 
cover in Colorado and there seem to be high-level clouds in the cirrus shield and 
frontal band casting shadows on the low-level clouds. When this happened, both the 
channel | and channel 3 reflectances were reduced considerably and the surface on 
which the shadow is cast was analyzed as land. This had a verv noticeable effect on 
the quatity of both of these analyses. 

The major problem with Case 3 was the analysis of cirrus. Three of the 
stations that were different from the surface observation were 4LJ (Lamar, CO), GCK 
Peaneemeciy 5), and LC (Hill City, KS), all due to the problemmvith cirrus. The 
discrepancy at LBF (North Platte, NE) appears to be a problem with the observation, 
perhaps an error in transnussion. The clouds in this area were handled verv well by the 
analysis and are to the southeast and west of LBF. Case 4 had discrepancies at 7 
stations all due to shadows on low clouds being cast by high clouds. Those stations 
are 4LJ (Lamar, CO), PUB (Pueblo, CO), TAD (Trinidad, CO), TCC (Tucumcari, 
NM), DHT (Dalhart, TX), GLD (Goodland, KS), and GCK (Garden City, KS). The 


69 







TABLE YS 


Analyzed versus reported sky cover for Case 3. 















Reported Analyzed Reported 
Station Coverage-21Z Coverage Coverage-22Z 
ALJ Scattered Clear Not available 
AKO Crean Clear Clear 
DDC Broken Broken Not available 
GCK Not Available Scattered Overcast 
GED Scattered Scattered Scattered 
ee Broken Clear Scattered 
DB Scattered Clear Overcast 
rex Overcast Scattered Scattered 





Clear 
Scattered 


Clear 


Broken 





Ee 
MCK 






Clear 


Overcast 









shadows were quite prominent in the visible image (see Appendix C, Fig. C.25) and 
also in the analysis shown here. The shadows were so prominent that the darker area 
in the analysis of snow cover between PUB and COS (Colorado Springs, CO) had the 
exact shape of the cloud element between PUB and LHX (La Junta, CO). Two other 
stations on Case 4 had discrepancies: COS and APA (Arapahoe, CO). For COS, low- 
level overcast was reported and the analysis was clear. The channel 4 brightness 
temperatures were mostly within 3 to 4 degrees of the reported surface temperature. 
There were some temperatures that were 10 degrees colder but these were in the higher 
elevations and the channel 3 reflectance associated with these pixels was 0.01. These 
characteristics are very indicative of snow cover. It is possible that there was a 
transmission or a data storage problem associated with this observation. The observed 
sky cover for APA was clear and the analysis was broken. The probability of an error 
in the observation is more likely here because there are quite obviously clouds in the 
visible and infrared images (see Appendix C, Figs. C.21 and C.23). The clouds that are 
present are even casting shadows so there may even be an overcast sky condition 


present at this station. Another explanation may hie in the fact that the 1500 GMT 
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TABLE 14 


Verification Case 3 cloud cover analysis results. 


Category Count Percentage 
Group 1] Correct 4/6 67 
l-category difference 1/6 17 
2-category difference 1/6 16 
| 3-category difference 0/6 0 
| Group 2 Correct 0/1 0 
l-category difference Val 100 
2-category difference 0/1 0 
3-category difference 0/1 0 
Group 3 Comect 2/3 6/7 
l-categorv difference 0/3 
2-categorv difference Oe 
3-category difference es 35 
Totals Correct 6/10 60 
l-category difference 2/10 20 
2-category difference 1/10 10 
3-categorv difference 1/10 10 


observation was not available for APA. This could be a case of clearing skies within 
the hour with the analysis being in between the changes. 
3. January 23, 1987 

The final two cases to be evaluated both came from NOAA-10 at 1443 GMT. 
A series of winter storms left a long swath of snow cover from New Mexico and Texas 
northeastward to the Great Lakes. Both images have snow-free land along with snow 
cover ranging from less than an inch in Iowa and Kansas to 11 inches in Indiana. The 
cloud cover in these images 1s only low- to mid-level clouds, with lake-effect clouds and 
snow showers in Michigan and Indiana. A complete sequence of channels 1, 3, and 4 
unages for these cases along with the analyses is shown in Appendix C, Figs. 
C.29-C.32 (Case 5) and C.33-C.36 (Case 6). 


7) 








Fig. 4.4 Same as Fig. 4.1 for Case 4. 


Case 5 shows nearly all of Kansas, the northern half of Oklahoma and 
northwestern Missouri. Of particular interest in this image is the low clouds directly 
over the snow cover in Kansas. Table 17 lists the observed sky cover for 1400 GMT 
and 1500 GMT, along with the analvzed sky cover. The results of the analysis are 
shown in Fig. 4.5 and the results of the verification are summarized in Table 18. Case 
6 covers eastern Iowa, southern Wisconsin and nearly all of Illinois, Indiana, Michigan 
and Lake Michigan. This image is from near the edge of the pass and has negligible 
snow cover except in Illinois and Indiana. This was a very cold airmass and the clouds 
in the image were generated by this very cold air flowing over the warmer Lake 
Michigan. Table 19 lists the observed sky cover for 1400 GMT and 1500 GMT, along 
with the analvzed sky cover. The analysis is shown in Fig. 4.6 and the verification 


results are summarized in Table 20. 


Case 5 1s perhaps the best of these 6 cases for illustrating the strength of this 
routine in handling low clouds over snow. The clouds in southern and western Kansas 
are apparent on the visible image only because of shadows around the edge of the 
cloud (Appendix C, Fig. C.29). They are virtually invisible on the infrared image 
‘Appendix C, Fig. C.31), except that they appear to be slightly darker than the cloud- 
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TAREE 1S 


Analyzed versus reported sky cover for Case 4. 


| Reported Analyzed Reported 
Station Coverage-15Z Coverage Coverage-16Z 
4LJ Broken Scattered Not available 
ALS Overcast Scattered Scattered 
AMA Broken Overcast Overcast 
APA Not available Broken Clear 
COS Overcast Clear Overcast 
SERS Overcast Clear Clear 
DHT Overcast Broken Overcast 
GCK Overcast Broken Overcast 
GLD Broken Scattered Overcast 
Ei Overcast Scattered Scattered 
ie Not available Clear Clear 
PUB Overcast Scattered Overcast 
SAF Scattered Broken Broken 
TAD Overcast Broken Overcast 
eC Overcast Broken Overcast 


free areas indicating the presence of a temperature inversion. The objective analysis 
cleared up the confusion quite dramatically and it even picked up a thin strip of cirrus 
cloud in Oklahoma that 1s evident on the infrared image but not on the visible image. 
Case 6 brought out the fact that the anisotropic factor is very critical when 
analvzing an image close to the edge of the pass. The only clouds in this image are 
directly over Lake Michigan with some cloud lines extending into Michigan and 
Indiana. These clouds can be seen fairly well in the images of Channels |, 3, and 4 in 
Appendix C, Figs. C.33 through C.35. The analysis overestimated the cloud cover, and 
It appears that this occurred over either clear land or land with very little snow cover. 
This is a morning pass so at the nght hand side of this image, the satellite zenith angle 


is large and the relative azimuth angle is large, along with a large solar zenith angle. 


73 


TABLE 16 


Verification of Case 4 cloud cover analysis results. 
| 
| 


Category Count Percentage 
Group | Correct 1/9 1] 
1-category difference oye 56 
2-category difference 2/9 Mp: 
3-category difference 1/9 1] 
Group 2 Correct 215 67 
l-category difference 1/3 33 
2-category difference 0/3 0 
3-category difference 0/3 0 
Group 3 Correct a8 100 
l-category difference 0/3 
2-category difference 0/3 
3-category difference 0/3 
Totals Correct 6/15 40 
l-category difference 6/15 40 
2-category difference pals 13 
3-category difference 1/15 ii 


At IND (Indianapolis, IN) the skies were reported clear and there was | inch of snow 
on the ground. The solar zenith, satellite zenith, and relative azimuth angles were 
about 74°, 60°, and 140° degrees respectively making the anisotropic factor 1.09. This 
factor is the average for bin 31 of snow and low cloud from the tables of Taylor and 
Stowe (1984) since land values were not reported for a solar zenith angle of greater 
than 72°. However, for a solar zenith angle ranging from 66° to 72°, the factor for 
land in this bin is 1.44. This value from the tables is based on only 6 observations, but 
is still much greater than the average that is used in the analvsis routine. What has 
happened is that the anisotropic factor used in the program was not as large as it 
should have been for this configuration of sun-satellite geometry. This caused the 


reflectance to be above the established threshold for land and resulted in many of these 
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Fig. 4.5 Same as Fig. 4.1 for Case 5. 


pixels being analvzed as cloud (see Eq. 2.12). Table 2.1 of Taylor and Stowe shows 
that land is more variable in terms of anisotropy than is snow, and values for land are 
not reported at a solar zenith of greater than 72°. It would have been better to have 
assumed that the values valid for the solar zenith angle range from 66° to 72° were also 
valid for the higher ranges for which values were not reported, rather than to just 
average the values in these ranges for snow and cloud. This would have produced a 
higher average value to be used by the program and mav have avoided some of this 
problem. 

The verification of Case 5 was disappointing because the analvsis appeared to 
produce excellent results. There was a discrepancy with all the stations in Oklahoma: 
ieieuisa), Tih {Tinker APB) END (Vance AFB) and PNC (Ponca Citv). The 
visible image shows surface features in the snow cover for END and PNC, and clear 
land is present for TIK and TUL (Appendix C, Fig. C.29). Channel 4 temperatures for 
all of these stations were within 3 to 5 degrees of the reported temperature, indicating 
no cloud cover at the height it was reported. This also appears to be the case at FOE 
(Forbes AFB, KS), where surface features are present in the snow cover on the visible 


image and were captured well by the analysis. At SLN (Salina, KS) and MKC (Kansas 
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TAD 


Analyzed versus reported sky cover for Case 5. 


Reported Analyzed Reported 
Station Coverage-14Z Coverage Coverage-15Z 
CNK Clear Clear Clear 
CNU Scattered Clear Clear 
EMP Not available Clear Clear 
END Overcast Clear Overcast 
POE Overcast Clear Overcast 
OTL Overcast Scattered Scattered 
JUN Overcast Clear leat. 
MHK Clear Clear Clear 
MKC Broken Overcast Broken 


PNC 
Roe 
5 1aN 
DER 

Tike 


Scattered 


Broken 
Clear 
Overcast 


Scattered 


Clear 
Broken 
Scattered 
Clear 
Clear 


Scattered 
Broken 
Cleat 

Scattered 

Scattered 





City, MO) the discrepancy is probably due to the different perspectives of the observer 
and satellite. At SLN the analysis had clouds only'‘at the outer limit of the sampled 
area making it very reasonable to assume the observer could not see them. At MKC it 
is reasonable to assume that the observer reported broken low-level clouds and was 
unable to see anything above these breaks. This overcast layer is evident in all images 
(see Appendix C, Fig. C.29-C.31) and was analyzed nearly perfectly by the routine. On 
Case 6, of the 15 stations with a discrepancy between the analysis and observation, 11 
can be attributed to the problem with the anisotropic reflectance factor. The other 4 
stations are GRR (Grand Rapids, MI), JXN (Jackson, MI), MKG (Muskegon, M1), 
and SBN (South Bend, IN). These stations reported greater cloud cover that was 
analyzed, possiblv due to clouds being directly over the station and the sampled area 


for each station including some clear areas. A 15-pixel radius was sampled in this 
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TAB EIESIS 


Verification of Case 5 cloud cover analysis results. 





Category Count Percentage 














Group | Sormect 4/10 40 
I-category difference 4/10 40 
2-category difference 0/10 0 
3-category difference 2/10 20 

Group 2 SOmrect 1/1 100 
I-categorv difference O/1 
2-category difference 0/1 
3-category difference 0/1 

Group 3 Correct 2/3 67 
I-category difference 1/3 We 
2-categorv difference 0/5 0 
3-category difference /3 0 

Totals SoOnicer 7/14 50 
l-categorv difference 5/14 36 
2-category difference 0/14 0 


3-category difference 2/14 14 






image and close to the right hand side, where all these stations are located, the 
resolution degrades to the point where this may be too large an area to be sampled. 
The discrepancy at MKE (Milwaukee, W1) can be attributed somewhat to the error in 
estimating the anisotropic factor as discussed above, but this is more likely due to most 
of the analyzed clouds being over Lake Michigan while directly over the station it was 
clear. 
4. Evaluation of the Analysis Routine 

The combined results of all six analyses are summarized in Table 21. It 1s 
especially encouraging that the results are heavily weighted on the positive side, that 1s, 
there were very few stations that had an analyzed sky cover two or three categories 
different from the analysis. 


The significant problems noted are as follows: 
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Fig. 4.6 Same as Fig. 4.1 for Case 6. 


Snow cover and cirrus have very simular reflective properties in channels | and 3 
so their correct analysis cannot be based on reflectance. There is a difference in 
the transmussive properties of cirrus between channels 3 and 4 that can be used 
to aid in the analvsis of cirrus, but there are still significant problems in that 
regard. It is estimated that about 50 percent of the cirrus on all of the images 
analvzed here was handled correctly based strictly on the channel 3 reflectance, 
with less than 5 percent added to this based on the channel 3-channel 4 
temperature factor. About 55 percent of allvcirrusiwas aidivzed canrect’* saimde 
temperature factor did not add very much to the analysis because solar 
reflection 1s included in the brightness temperature difference between channels 
3 and 4. The attempt was to analyze cirrus based on a brightness teinperature 
difference due the the different transmissive properties of “cirrus between 
channels 3 and 4, but the effects of solar reflection nearlv completely wiped out 
any detectable differences in thermal emission. Cirrus clouds only within a 
narrow range of optical depths could be unambiguously separated from snow 
Coy ae 

Channel | reflectance was used in the analysis of land. There were instances 
where there was some overlap with hquid clouds and the threshold setting was 
such that some land was analyzed as cloud. Much more significant in the 
analysis of land was the fact that anisotropic factors are not available for land 
when the solar zenith angle is greater than 72°. This was critical because 
according to the values of Taylor and Stowe (1984); land is more variable in 
terms of anisotropy than is snow or cloud. In Case 6, clouds were incorrectly 
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Station 


ALO 
AZO 
BMI 
BRL 
CID 
CMI 
CWA 
DBO 
DEC 
DSM 
EAU 
FWA 
GRB 
GRR 
GUS 
HUF 
IND 
IXN 
LAF 
LAN 
LNR 
LSE 
MCW 
MKE 
MKG 
MLI 
MSN 


TABLES 


Analyzed versus reported sky cover for Case 6. 


Reported Analyzed Reported 
Coverage-14Z Coverage Coverage-15Z 
Clear Clear Clear 
Scattered Broken Scattered 
Clear Clear Clear 
Clear Clear Not available 
Clear Clear Not available 
Clear Clear Not available 
Clear Glcar Not available 
Clear Clear Clear 
Clear Clear Clear 
Clear Clear Clear 
Clear Clear Clear 
Scattered Broken Not available 
Clear Scattered Clear 
Overcast Scattered Overcast 
Overcast Broken Glear, 
Clear Clear Clear 
Clear Scattered Clear 
Overcast Scattered Overcast 
Clear Scattered Clear 
Overcast Scattered Clear 
Clear Clear Not available 
Clear Scattered Clear 
lear Clear Clear 
Clear Broken Clear 
Overcast Broken Overcast 
Clear Clear Clear 
Clear Scattered Clear 


TABLE 719 


Analyzed versus reported sky cover for Case 6. (cont’d.) 


Not available Broken Overcast 
Clear Scattered Clear 
Clear Clear Clear 
Clear Clea Clear 
Clear Scattered Clear 
Clear Clean Clear 

Overcast Broken Overcast 


Not available Clear Clear 





analyzed over clear land close to the edge of the image because the anisotropic 
factor used was less than what it should have been. 


An unexpected problem was with shadows. The cloud analysis was affected by 
this when high-level clouds cast shadows on low-level clouds. For the reader 
who 1s more interested in a correct analysis of land versus snow, this problem 
affected the analysis of snow cover. The shadows caused the reflectance in both 
channels 1 and 3 to be reduced and the affected pixels were analyzed as land. 
This problem was particularly noticeable in the winter images when the solar 
zenith angle was the largest. Unlike the other problems, nothing can be 
objectively done about this. 


There are also some reasons for discrepancies between the analyzed and 


observed sky cover that are beyond the control of the routine: 


ie 


The time difference between the surface observation and the satellite pass was 
great enough in some cases to lead to a significant difference in sky cover. This 
occurred when the weather system was very dynamic and there were rapid 
changes in cloud cover. It is for this reason that the analvses were compared 
with the hourly observations on either side. 


The perspectives of a ground observer and a satellite are very different, the 
ground observer viewing the clouds from below and the satellite viewing the 
clouds from above. The ground observer does not spend hours meticulously 
evaluating the exact coverage of each individual cloud element and then add 
these together to arrive at a total cloud coverage. A computer can take the 
time to do this but the observer makes a subjective determination of both 
height and sky cover and no two observers will make the same determination at 
anv given time. In a convective situation much of the reported sky cover is due 
to the sides of the clouds, especially those viewed at a large zenith angle. 
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TABLE 


Verification of Case 6 cloud cover analysis results. 





Category Count Percentage 


Group } Correct 18/33 55 
l-category difference PB 8) 36 
2-category difference eo 9 
3-category difference Oras 0 
Group 2 Correct 0/0 0 
]-category difference 0/0 0 
2-category difference 0/0 0 
3-category difference 0/0 0 
Group 3 Gorect 2/2 100 
]-category difference 0/2 0 
2-category difference 0/2 
3-category difference 0/2 0 
Totals Correct 20/35 57 
]-category difference 1235 34 
2-category difference 31D S 
3-category difference 0/35 0 





Because of this, gaps between the cloud elements become smaller and disappear 
as the zenith angle increases. The resulting sky cover will be greater than what 
a satellite would see for that station since all the clouds within the sampled 
radius can be viewed within just a few degrees of scan angle. An observer may 
be limited to the portions of the skv that can be evaluated due to obstructions 
around the airfield and, in addition, an observer may intentionally limit the 
distance from the airfield at which cloud cover will be evaluated, concentrating 
onlv on what 1s in the immediate vicinity of the airfield. 


tod 


Characteristic to all polar-orbiting satellites 1s errors in navigation. Error in 
location of the individual stations on the images was about 3 to 5 kilometers. 
The geographical boundaries were corrected on images where landmarks were 
visible for assistance, but no attempt was made to correct the location of the 
individual stations. 
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TAHIBE 21 


Combined statistics for all analyses. 


Category Count Percentage 
Carect 60/110 58 


l-category difference 36/110 33 


2-category difference 10/110 9 


3-category difference 4/110 





4. Errors in the surface observations arose due to problems with storage in the 
data base or errors in transmission. Sky cover in the observations used in this 
study is stored in a single character, making the probability of an error greater 
than if it were stored as the standard three-character abbreviation. Many 
stations had to be discarded as potential candidates for verifying the analyses 
because meaningless characters for sky cover had crept into the data base. 


It is unfortunate that so many of the discrepancies between the analysis and 
observations were attributed to problems beyond the control of the routine. It would 
have been much more tedious to get copies of actual coded observations from the 
archives, but that most certainly would have yielded improved verification results. In 
spite of these problems, however, the statistics indicate that the analysis routine was 


very successful under a wide variety of circumstances. 
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V. SUMMARY AND RECOMMENDATIONS 


This thesis reported values of solar reflectance of snow, clouds and land at 3.7m 
using data from channel 3 of the AVHRR. The reflectance was derived using a 
method of removing the thermal emission portion of the channel 3 radiance 
measurement. These values were calculated using a rough estimation of anisotropic 
reflectance and values agree quite well with those determined from various theoretical 
studies. The reason is that for most clouds, thermal emission is a much smaller part of 
the channel 3 radiance measurement than js solar reflection (see Fig. 2.3), so any errors 
in the technique that removes thermal emission introduce little error into the derived 
reflectance. The greatest error was noted for the reflectance of snow cover. Its 
reflectance is typically less than 5 percent and in this situation solar reflection and 
thermal emission are nearly equal contributors to the channel 3 measured radiance. 
Thus the error in removing the thermal emission translated into larger error in 
reflectance. 

An automated cloud analysis was then developed to handle the significant 
problem of snow/cloud discrimination. The routine exploits the fact that snow has a 
very low reflectance in AVHRR channel 3 while liquid clouds and most ice clouds have 
a high reflectance. Since the reflective properties of snow cover and many ice clouds 
are sO similar in both channels 1 and 3, the routine attempted analysis of cirrus that 
did not meet the channel 3 reflectance threshold based on a difference in its 
transmissive properties between channels 3 and 4. Unfortunately, this technique did 
not work well because in addition to the difference in transmissive properties, solar 
reflection contributes to the difference in brightness temperature. The reflectance was 
variable enough to make unambiguous separation of the two surfaces possible onlv for 
cirrus of a limited range of optical depths. The algorithm was very successful and 
produced excellent results, even at very large solar zenith angles. Its main strength was 
in separating low clouds from snow cover, while its weakest point was the analvsis of 
cirrus clouds. The analyses were verified against the coincident surface observations 
and 88 percent of the analyzed stations either matched or were just one category of sky 


cover different from the observed sky cover. 
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Although the routine produced generally excellent results, some improvements 


are recommended: 


ie 


tJ 


It was suggested in Chapter I] that additional theoretical studies be undertaken 
to improve the available information on the anisotropic reflectance factor. In 
the interim, what we know can be put to better use. First needed 1s a verv easy 
near-term solution to the missing values for land at solar zenith angles greater 
than 72°. The factor currently used is the average for cloud and snow. The 
land values for the range from 66° to 72° could be used for the higher ranges, 
thus vielding a higher average for use by the analysis routine. Since the 
configuration of NOAA satellites limits the range of satellite zenith angle and 
thus the number of bins that need to be used by the computer program, more 
accurate values can be used as input to the program. Rather than the 
averaging technique used in this study, a table of values for each surface should 
be made available for the program to choose from. A preliminary guess of the 
scene classification could be made and then an iterative process used to arrive 
at the correct factor for that particular classification. 


For the most part, cirrus could be separated successfully by this routine by 
setting the channel 3 threshold at 0.057. Approximately 50 percent of the cirrus 
was analyzed correctly with this threshold. An attempt was made to improve 
the analysis of the cirrus that was too thin to be separated with this threshold 
by looking at the difference in the transmissive properties between channels 3 
and 4. In the presence of cirrus, there is a difference in brightness temperature 
between these two channels while there ts little difference in the clear sky case. 
This was formulated into a unitless factor and then a threshold was developed 
from the test data. This factor did not add significantly to the analysis of cirrus 
because variable solar reflection “contaminated” the difference in brightness 
temperature due to thermal emission. This characteristic of cirrus in these two 
channels should be explored further to see if any additional techniques can be 
developed to improve the analysis. Another possibility is to work out a way to 
incorporate surface temperature into the analysis routine. A statistical analysis 
of all the clear sky pixels in the image could be used to determine a dynamic 
threshold temperature based on conditions present in that particular image. 
Development of a technique of this kind will require careful evaluation of 
computer processing time. As it is now written, the routine requires 20 to 30 
minutes to generate a single analysis so it is already computationally-intensive. 
In spite of added quality, the additional time required may lessen the usefulness 
of the analysis. 


More work should be done to compile a data base of near-infrared reflectances. 
The values reported here were based on samples of 4800 pixels from two 
satellite passes. This may seem to be a large number of pixels but it is not an 
extensive sampling. The sampling was from only two days of data in one 
season of the vear and with a very limited range of sun-satellite geometry. 
Building a complete data base should span several future studies and work must 
continue so that software can be developed to process data from the 1.6m 
sensors when thev become operational. 
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This analysis routine is a medium-sized software package written in FORTRAN 
language and it is adaptable to any computer system. However it is a computationally- 
intensive routine. Each of the analyses took from 20 to 30 minutes to generate on a 
Digital Equipment Corporation VAX 11/780 computer system. This compared with 2 
to 5 minutes needed to generate an image from any of the individual channels. It 
should be pointed out that some of this time was due to the program reading the 
satellite data from tape, but there are many calculations done to derive the channel 3 
reflectance. This routine could be put to use directly on a minicomputer system with a 
graphics capability and would generate analyses in a reasonable time period for 
operational purposes provided enough disk storage were available for the data, and 
provided that the scale were reduced from what was presented here. RTNEPH cannot 
use this routine directly as written because the processing is so much different. 
RTNEPH processes the visible satellite data and makes a cloud decision based on a 
comparison of reflectance at an analysis point with a climatological background 
reflectance for that point. If the analysis has a greater reflectance than the background 
value, a positive decision 1s made. Of course this procedure fails when the background 
1s snow because there is no difference in reflectance. Adding a component that 1s 
called only for points that are snow-covered could be the solution to this problem. 
This component would process near-infrared reflectance from the satellite data and 
compare this with a background near-infrared reflectance. This study confirms the 
verv low near-infrared reflectance of snow compared with the relatively high near- 
infrared reflectance of clouds so this approach could be used with a great deal of 
success, however, implementation of this depends on the availability of a data base of 


near-infrared reflectances as discussed above. 
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APPENDIX A 
SYMBOLS AND CONSTANTS 


SYMBOLS 
Ma SO a So ak cw os oe solar zenith angle 
Oc eee whales hace en nn er ee satellite zenith angle 
(i ee ee So ca na oot relative azimuth angle 
Oe OM so ee cad observed radiance 
BCT) . os isle en ve 4 eRe G4 ee es ce Planck function 
| are ra OR i opp cen soo & incident solar radiance 
er i ks 0 solar spectral irradiance 
I i EE ci EC FG oS CU Uae Oe eb 8 emissiVitV 
ee ee cc Oe 2 OM cone reflectance 
1(0 0,0) ee. 5 25 RMI co o clsleee clet directional reflectance 
er ere er eM so oo eS transmissiVitv 
re nr nt. of ee a so hee 6.0 0 ay anisotropic reflectance factor 
WR oO od 5 ba gee d 3 percentage of a pixel covered by clear land 
Wen oo een Otis a pau a percentage of a pixel covered by liquid clouds 
Wc eee Goh cab ag) wedge percentage of a pixel covered by ice clouds 
WME nr a came ene cb och esas anc percentage of a pixel covered by snow 
| ee cob < acneyielin o oo © wavelength 
V oo s Gulla He ee © He aoe One gee chunk Oe nn wavenumber 
To aw wile BG ale paw ties fo dace ate hae eter ete Pree brightness temperature 
Fo eee oe eee channel 3/channel 4 brightness temperature factor 
CONSTANTS 
h = Planck s Comstantis 004... 20 ae een nee ere 6.6262 x 10°74 Joule sec 
K = Boltzmann s Constant... ae eee nee eee escOomecmlOn. Joule/deg 
c = velocity of light--4075 <a eee << pee oe 10° meter, sec 
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APPENDIX B 
CALCULATION OF AVHRR CHANNEL 3 RADIANCE 


Fig. B.1 shows the solar spectral irradiance for the channel 3 spectral band, from 
measurements taken by Thekaekara er al. (1969) at an altitude of 11.58 kilometers. 
Irradiance values were converted to wavenumber units since channel 3 AVHRR data 
are calibrated in those units. Since there is no significant absorption below that 
altitude, these values were assumed to be the incident solar irradiance at the earth's 
surface. Linear interpolation was used to give values at wavenumbers for which 
response functions are reported. 

The following equation was used to calculate the normalized irradiance for 


channel! 3: 


N N 
E,(v) = ¥ E(v,)9(v,Av,/> @(v, Av, (B.1) 


i= 1 i=] 


where @(V.) is the spectral response for each interval in the band. Figs. B.2 and B.3 
show the normalized spectral response for NOAA-9 and NOAA-10 channel 3 taken 
from Lauritson es al. (1979). Assuming isotropic radiation, the solar radiance was 


found by: 
[,(v) = E,(v)/n (Bae) 


Using these equations the normalized radiance for NOAA-9 is 5.31085 
mW m?.cm’].sr, and for NOAA-I0O the normalized radiance is 5.26415 mW/ 


A =) 
MmMeCtn «ST. 
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Fig. B.} Solar spectral irradiance for NOAA-9 and NOAA-10 channel 3 spectral bands. 
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Channel 3 (NOAA-9) normalized spectral response. 
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Fig. B.3. Channel 3 (NOAA-10) normalized spectral response. 
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APPENDIX C 
SATELLITE IMAGES 


This appendix combines all the images that were analyzed and verified in this 
thesis in one convenient location. Included with all analyses are the coincident visible 
and infrared images for comparison. Also included is an image of channel 3 brightness 
temperature. The analyses that were presented in Chapter IV are presented again here 
and the snow-covered areas are enhanced to show the difference between snow-covered 
and clear land. This illustrates the advantage of displaying land, snow, and clouds on a 
unique range of grayshade values. Any range of values can be enhanced as desired by 
the analyst and even various color enhancement schemes can be easily employed by 


displaying in this manner. 
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Fig. C.l 9 Nov 86 subscene 1, NOAA-9 AVHRR channel | reflectance. 





Fig. C.2 9 Nov 86 subscene 1, NOAA-9 AVHRR channel 3 brightness temperature. 
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Fig. C.4 9 Nov 86 subscene 1, results of cloud cover analysis. 
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Fig. C.8 Same as Fig. C.4 for 9 Nov 86 subscene 2. 
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Fig. C.12 Same as Fig. C.4 for 10 Nov 86 subscene 1. 
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Fig. C.21 Same as Fig. C.1 for Case 3. 





Fig: C.22 Samerasehigs Geto, Cases. 
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Fig. C.24 Same as Fig. C.4 for Case 3. 
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Fig. C.26 Same as Fig. C.2 for Case 4. 
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Fig. C.28 Same as Fig. C.4 for Case 4. 


105 


Fig. C.30 Same as Fig. C.2 for Case-5. 
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Fig. C.34 Same as Fig. C.2 for Case 6. 
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Same as Fig. C.3 for Case 5. 
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